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ABSTRACT 

In a dynamical-radiative model we recently developed to describe the physics of compact, GHz- 
Peaked-Spectrum (GPS) sources, the relativistic jets propagate across the inner, kpc-sized region of 
the host galaxy, while the electron population of the expanding lobes evolves and emits synchrotron 
and inverse- Compton (IC) radiation. Interstellar-medium gas clouds engulfed by the expanding lobes, 
and photoionized by the active nucleus, are responsible for the radio spectral turnover through free- 
free absorption (FFA) of the synchrotron photons. The model provides a description of the evolution 
of the spectral energy distribution (SED) of GPS sources with their expansion, predicting signifi- 
cant and complex high-energy emission, from the X-ray to the 7-ray frequency domain. Here, we 
test this model with the broad-band SEDs of a sample of eleven X-ray emitting GPS galaxies with 
Compact-Symmetric-Object (CSO) morphology, and show that: (i) the shape of the radio continuum 
at frequencies lower than the spectral turnover is indeed well accounted for by the FFA mechanism; 
(ii) the observed X-ray spectra can be interpreted as non-thermal radiation produced via IC scatter- 
ing of the local radiation fields off the lobe particles, providing a viable alternative to the thermal, 
accretion-disk dominated scenario. We also show that the relation between the hydrogen column 
densities derived from the X-ray (Ah) and radio (Nm) data of the sources is suggestive of a positive 
correlation, which, if confirmed by future observations, would provide further support to our scenario 
of high-energy emitting lobes. 

Subject headings: galaxies: active - galaxies: individual (IERS B0026+346, IERS B0108+388, 
IERS B0500+019, IERS B0710+439, PKS B0941-080, IERS B1031+567, IERS 
B1345+125, IVS B1358+624, IERS B1404+286, IERS B2128-I-048, IERS 
B2352+495) galaxies: jets - X-rays: galaxies - radiation mechanisms: non- 
thermal 



1. INTRODUCTION 

The power released by active galactic nuclei (AGNs) is 
currently interpreted in terms of conversion of gravita- 
tional energy to radiative energy by accretion processes 
feeding the central supermassive black hole (BH) with 
environmental gas. The triggering, maintenance, and 
fading of the AGN activity, as well as the link of these 
processes with the physical conditions of the environ- 
ment, from sub- to super-galactic scales, are still widely 
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debated issues, and keep on stimulating a large variety 
of scientific investigations. In this context, radio galaxies 
are ideal laboratories, because they offer an edge-on view 
of both their nuclei and their relativistic jets, launched 
from the galactic center and reaching up to Mpc dis- 
tances. The variety of powers, sizes, and morphologies 
displayed by radio galaxies not only provides pieces of 
evidence of different environmental physical conditions, 
but also samples subsequent stages of the source evolu- 
tion. In particular, key sources for the investigation of 
the very first phases of the evolution of radio galaxies 
are the Gigahertz-Peaked-Spectrum (GPS) sources as- 
sociated with galaxies and characterized by a Compact 
Symmetric Object (CSO) radio morphology. 

GPS sources (see O'Dea 1998, for a review) are a class 
of powerful radio sources (Pi.4GHz ~ 10 25 W Hz -1 ) dis- 
playing convex radio spectra that turn over at frequencies 
of about 0.5-1 GHz; they make up a conspicuous fraction 
(~10%) of the radio sources found in high-frequency ra- 
dio surveys, and are optically identified with either galax- 
ies or quasars; their radio sizes are smaller than about 
1 kpc, and their radio morphologies reveal either core- 
jet structures or mini-lobes embedding terminal hotspots 
and possibly straddling a central core. CSOs (Wilkinson 
et al. 1994; Readhead et al. 1996) have instead emerged 
in VLBI surveys as a class of radio sources with compact 
(^500 pc) and symmetric radio structures, making them 
resemble "classical double" radio galaxies in miniature. 
Whereas there is no debate on the true compactness of 
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CSOs, the compactness of a fraction of the core-jet GPS 
sources might be the result of foreshortening of extended 
radio sources roughly aligned with the line of sight; espe- 
cially (although not exclusively) in these cases, the GPS 
itself might be a transient spectral state: flux-density 
variability and polarization studies at different radio fre- 
quencies are thus instrumental to the selection of bona- 
fide samples of GPS sources (Tinti et al. 2005; Torni- 
ainen et al. 2005, 2007; Orienti & Dallacasa 2008a). The 
overlap between the GPS-source and the CSO classes is 
however significant: GPS sources associated with galax- 
ies are most likely to display a CSO morphology, and 
most CSOs exhibit a GPS. This overlap is explained in 
terms of synchrotron radio spectra dominated by the 
emission of the mini-lobes, and suffering from absorp- 
tion effects causing the turnover about 1 GHz (Snellen 
et al. 2000) . GPS /CSOs associated with galaxies are thus 
high-confidence candidates for truly compact sources. 

Although compact GPS sources were proposed to be 
young objects soon after their discovery (e.g., Shklovsky 
1965), a widely discussed alternative to the youth sce- 
nario was the frustration scenario, ascribing the small 
source size to confinement by a particularly dense inter- 
stellar medium (ISM) preventing the jet expansion (van 
Breugel et al. 1984). Because the required ISM confin- 
ing densities are not confirmed by recent studies (e.g., 
Morganti 2008), the confinement scenario is considered 
less likely, although it might still apply to selected ob- 
jects (e.g., Garcia-Burillo et al. 2007). Much observa- 
tional evidence has instead accumulated in favour of the 
youth scenario, the most compelling measurement being 
the detection, in several GPS/CSO galaxies, of hotspot 
advance velocities about 0.1-0. 2c (Owsianik & Conway 
1998; Owsianik et al. 1998; Tschager et al. 2000; Taylor et 
al. 2000; Gugliucci et al. 2005), which indicate kinemat- 
ical source ages not higher than a few 10 3 years, in good 
agreement with spectral ageing estimates (Murgia et al. 
1999). Further evidence of youth comes from the under- 
luminosity of GPS/CSO optical narrow-emission lines, 
suggesting that the Stromgrcn sphere of the recently- 
triggered AGN is still in an expansion phase in these 
sources, and we are thus witnessing the birth of their 
Narrow-Line Region (NLR; Vink et al. 2006). 

The youth scenario is part of a wider evolutionary sce- 
nario (Fanti et al. 1995; Snellen et al. 2000), according 
to which GPS/CSOs would first evolve into symmetric 
Compact-Steep-Spectrum (CSS) sources, equally power- 
ful but with sizes of about 1-15 kpc, and then further 
expand outside the host galaxy to become large-scale, 
powerful radio galaxies. The evolutionary link between 
GPS and CSS sources is supported by their membership 
to the anticorrelation between linear size and peak fre- 
quency (LS oc i/p 0,65 : O'Dca & Baum 1997). However, 
it is not clear yet whether all GPS/CSOs will eventu- 
ally become large-scale radio sources, and whether the 
evolved sources will display a Fanaroff- Riley type I (FRI) 
or type II (FRII) morphology, because of the frequency's 
observational limits in the exploration of this relation for 
super-galactic sized sources, as well as because of the 
biases in the surveys from which the distribution of the 
radio sources in the power- linear size (P — LS) diagram is 
drawn (e.g., Snellen et al. 2000, and references therein). 

The first systematic studies of the host galaxies of 



GPS sources, conducted in the optical and near-infrared 
(NIR) bands (Snellen et al. 1996; O'Dea et al. 1996; de 
Vries et al. 1998a,b, 2000) showed that they are, as CSS 
sources, characterized by luminous masses, brightness 
profiles, and optical-NIR colors more typical of passively- 
or non-evolving giant elliptical galaxies, than either spi- 
rals, small ellipticals, brightest cluster galaxies (BCGs), 
or central-dominant (cD) galaxies at similar redshifts. 
This seems to be consistent with what was found for 
a sample of FRIIs rather than with the properties of 
FRIs, which are preferentially hosted by cD galaxies 
(Owen & Laing 1989; Owen & White 1991; Zirbel 1996), 
and suggest that GPS sources are more likely to evolve, 
through the CSS phase, in FRIIs rather than in FRIs (de 
Vries 2003). However, the morphologies of the majority 
(~60%; O'Dea et al. 1996) of the hosts of GPS sources 
show signs of recent mergers and/or interactions. This 
evidence, apparently at odds with the passively-evolving 
scenario, was first interpreted as an indication that GPS 
sources might be associated with the first of a sequence 
of mergers (de Vries 2003), possibly leading to the for- 
mation of a BCG or a cD only over time scales much 
longer than the lifetime of the radio source. In sub- 
sequent investigations, however, the presence of young 
stellar population in these giant ellipticals was revealed 
by means of stellar-population synthesis models applied 
to photometric (de Vries et al. 2007) and spectropho- 
tometric data (Holt 2009). Furthermore, de Vries et 
al. (2007) showed the similarity, out to z ~ 0.55, be- 
tween the luminosities of the GPS-source hosts and of 
the Luminous Red Galaxies (LRGs), which are thought 
to represent the most massive early-type galaxies, and 
are often associated with BCGs. Therefore, the proper- 
ties of the hosts do not provide unambiguous indications 
on the subsequent GPS/CSO evolutionary stages. 

Despite these results, population studies highlight the 
existence of far too many compact sources compared to 
the number of powerful large-scale objects (O'Dea & 
Baum 1997). Beside the possibility, not easy to justify, 
that these sources undergo a short-lived activity (Read- 
head et al. 1996), two scenarios prove to be particularly 
attractive to explain this statistical evidence. In the first 
scenario, the jet activity is intermittent on time scales 
comparable to the age of the small sources (e.g., Reynolds 
& Begelman 1997): this scenario, supported by the obser- 
vations of double-double radio galaxies (Lara et al. 1999; 
Schoenmakers et al. 2000; Kaiser et al. 2000), and by the 
detection of candidates for dying compact sources (Giro- 
letti et al. 2005; Parma et al. 2007), finds a natural ex- 
planation in the framework of accretion-disk instabilities 
(Czerny et al. 2009). An accretion disk operating above 
a given threshold accretion rate (m* ~ 0.025 rriEdd) is 
expected to experience instabilities driven by the ra- 
diation pressure, causing the source undergo alternate 
phases of high and low activity; for moderate accretion 
rates, the duration of the active phases is short enough 
(~ 10 3 — 10 4 years) to make the compact source unable to 
grow to super-galactic sizes, whereas accretion rates close 
to the Eddington limit are required for the development 
of large-scale sources (Czerny et al. 2009). In the sec- 
ond scenario, many sources undergo processes of jet-flow 
turbulent disruption before their lobes can grow to large 
sizes, either disappearing from the radio sky (Alexander 
2000) or developing an FRI-typc morphology (Kaiser & 
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Best 2007), in both cases experiencing a drop in lumi- 
nosity. The recent finding that FRIs smaller than ~40 
kpc were almost absent in a sample of radio galaxies op- 
tically identified in the SDSS (Best 2009) seems to sug- 
gest that all radio sources begin their life with collimated 
jets, but less powerful jets are more easily disrupted, as 
the source grows, in denser environments, leading to the 
production of FRI sources. Further evidence supporting 
this view might derive from the confirmation that the 
large fraction of low-power, compact objects with non- 
FRI morphology discovered within a sample of FRI can- 
didates at 1 < z < 2 (Chiaberge et al. 2009), are indeed 
young sources. The two above scenarios are not nec- 
essarily in conflict with each other: for sufficiently high 
accretion rates, accretion-disk instabilities might partake 
to the radio-source evolution, otherwise driven by the in- 
terplay between jet power and environmental conditions. 
However, this picture needs to be further explored. 

It has recently become clear that instrumental clues on 
the evolution of young GPS sources may come from the 
X-ray electromagnetic window. Early detections of GPS 
sources by ASCA (O'Dea et al. 2000) triggered extensive 
searches for X-ray emission from these sources; thanks 
to the capabilities of XMM-Newton and Chandra, several 
GPS/CSOs are now known to be strong X-ray emitters 
(Guainazzi et al. 2004; Vink et al. 2006; Guainazzi et al. 
2006; Siemiginowska et al. 2008; Tengstrand et al. 2009). 
However, the origin of this X-ray emission is not known 
yet. 

The best spatial resolution currently available in the 
X-ray band (~1" with Chandra) is not sufficient to re- 
solve the X-ray morphology of most GPS/CSOs: thus, 
the identification of the origin of the X-ray emission in 
these objects relies entirely on spectral studies. The com- 
ponents of GPS/CSOs first proposed to be the source 
of the observed X-rays are the ISM of the host galaxy 
shocked by the expanding radio lobes (Heinz et al. 1998; 
O'Dea et al. 2000), and the accretion disk's hot corona 
(Guainazzi et al. 2004; Vink et al. 2006; Guainazzi et al. 
2006; Siemiginowska et al. 2008): both components are 
expected to generate thermal X-ray radiation, but they 
are difficult to disentangle and do not rule out alterna- 
tive scenarios (Siemiginowska et al. 2009, and references 
therein). Furthermore, the location of GPS/CSOs on the 
radio/X-ray luminosity plane does not fully elucidate the 
relationship between "young" and "mature" radio galax- 
ies: GPS/CSOs are as X-ray luminous as FRIIs, despite 
their larger radio luminosity, but they seem to lie on the 
high radio-power tail of the radio-core/X-ray correlation 
discovered for FRIs (Tengstrand et al. 2009, and refer- 
ences therein). 

A way of investigating the origin of the X-rays in 
GPS/CSOs, and address, at the same time, some of 
the open issues on the physics and fate of young radio 
sources, is to simultaneously model the dynamical and 
radiative evolution of GPS/CSO galaxies, and constrain 
the models through the wealth of multiwavelength data 
currently available. We recently proposed a dynamical- 
radiative model that, for the first time, describes the evo- 
lution of the broad-band emission of a young GPS with 
a given jet kinetic power, as it expands through the ISM 
of the host galaxy (Stawarz et al. 2008). The model ac- 
counts for the radiative contribution of the various AGN 
components, as well as for environmental absorption ef- 



fects, and predicts significant and complex non-thermal 
X-ray to 7-ray emission. 

In the present paper, we apply the model to a sample of 
GPS/CSO galaxies, with the manifold aim of reproduc- 
ing their complex radio to X-ray spectra, and constrain- 
ing relevant source parameters as the jet kinetic power, 
the accretion rate, and the absorption mechanisms. New 
observational evidence supporting the model is also dis- 
cussed. The paper is organized as follows: in Section 
2, we briefly review the main features of the dynamical- 
radiative model we presented in Stawarz et al. (2008); in 
Section 3, we describe the source sample that we chose 
for the application of the model, and the relevant broad- 
band data; in Section 4, we show the results of the SED 
modeling; in Section 5, we present further evidence sup- 
porting the proposed scenario; we discuss our results in 
Section 6, and in Section 7 we draw our conclusions. 

A model of A-dominated universe (J1a = 0.7, Qm = 
0.3, il k = 0; Spergel et al. 2003), with H Q = 72 
km s _1 Mpc -1 (Freedman et al. 2001) will be adopted 
throughout this paper. Source-intrinsic quantities from 
the literature given below were corrected for this cosmol- 
ogy, unless otherwise stated. 

2. THE MODEL 

We recall below the most relevant features of our 
non self-similar dynamical-radiative model, presented in 
Stawarz et al. (2008). The reader is referred to the orig- 
inal paper (and references therein) for a more compre- 
hensive discussion. 

2.1. Dynamical Evolution 

We based our description of the dynamical evolution 
of a young radio source in its "GPS phase" on the model 
proposed by Begelman & Cioffi (1989) to explain the ex- 
pansion of classical double sources in an ambient medium 
with density profile p(r). The relevant equations were de- 
rived by assuming that: (i) the jet momentum flux, which 
is proportional to the jet kinetic power Lj, is balanced 
by the ram-pressure of the ambient medium spread over 
an area A^; (ii) the lobe sideways expansion velocity, v c , 
equals the speed of the shock driven by the overpres- 
sured cocoon, with internal pressure p, in the surround- 
ing medium; (iii) the kinetic energy transported by the 
jet pair during the entire source lifetime t is converted 
into the cocoon's energy pV (where V is the volume of 
the cocoon) at the heads of the jets. 

For a young GPS source expanding in the central 
gaseous core of a giant elliptical galaxy, and character- 
ized by age t, linear size LS(t) ^ 1 kpc, and transverse 
size l c (t), we could constrain the model with a number of 
reasonable approximations: (i) a constant ambient den- 
sity p = m p n (with m p the proton mass, and n — 0.1 
cm~ 3 ), representative of the inner King density profile 
of the host galaxy; (ii) a constant hot-spot advance velo- 
city Vh, as suggested by many obervations of CSOs (see 
Section 1); (iii) a scaling law l c (t) ~ i 1 / 2 , reproducing 
the initial, ballistic phase of the jet propagation, accord- 
ing to Kawakatu & Kino (2006) and Scheck et al. (2002) 
(see, however, Kawakatu et al. 2009a, for an alternative 
scenario). Under the aforementioned conditions, all the 
lobe physical quantities become functions of two param- 
eters only: the jet kinetic power Lj , and the source linear 
size LS. 
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2.2. Spectral Evolution 

In the framework of the dynamical model described 
above, we studied the evolution of the broad-band radia- 
tive output of GPS sources with the source expansion, for 
a given jet kinetic power Lj . The magnetic field in the ex- 
panding lobes scales as B = (Stttibp) 1 ^ 2 ~ L J 1//4 L5' _1/ ' 2 , 
with ?7b = Ub/p < 3, and Ub being the magnetic en- 
ergy density. The electron population (5(7) (with 7 the 
electron Lorentz factor), injected from the terminal jet 
shock to the expanding lobes, evolves under the joint ac- 
tion of adiabatic and radiative energy losses, yielding a 
lobe electron population N e (jy), characterized by a bro- 
ken power-law form with critical energy 7 cr when Q(7) 
is a power law, and by a more complex form when (5(7) 
is a broken power law with intrinsic break 7; nt ~2x 10 3 
(7int = mp/me, with m p and m e the proton and elec- 
tron masses, respectively; Stawarz et al. 2007). Assuming 
that the lobe electrons, in rough cquipartition with the 
magnetic field and the protons (e.g., Orienti & Dallacasa 
2008b, and references therein), provide the bulk of the 
lobe pressure, the electron energy density is U c = r/ c p, 
with i] c ^ 3. 

The lobe electrons emit synchrotron radiation, with 
luminosity L syn roughly constant with time, and en- 
ergy density U syn ~ LS~ 3 ^ 2 . This radiation suffers then 
from absorption processes, responsible for the character- 
istic spectral change across the turnover at GHz frequen- 
cies. Because the coupling between the synchrotron-self- 
absorption (SSA) process and our model equations re- 
turns turnover frequencies v v systematically lower than 
those typically observed, and reproduces neither the 
observed slopes of the optically-thick spectra, nor the 
v p — LS anticorrelation, the absorption process that we 
favoured and implemented is the free-free absorption 
(FFA). 

Recent data on compact sources of different brightness 
appear to be in good agreement with the predictions of 
the SSA model for homogeneous sources (de Vries et al. 
2009). However, this result does not exclude the viability 
of the FFA models (Bicknell et al. 1997; Begelman 1999). 
In fact, de Vries et al. (2009) showed that the Bicknell et 
al.'s FFA model works well on bright sources, although 
it does not on faint sources, unless the ratio between 
jet power and radio luminosity depends on luminosity. 
Although this scenario might be more complicated than 
homogeneous SSA, this solution was never investigated 
either in Bicknell et al.'s or in Begclman's FFA models 
in sufficient detail. In any case, this issue is not of ma- 
jor relevance for our sample, because we focus on bright 
sources, for which the FFA model works properly. 

In particular, we adopted the FFA scenario by Begel- 
man (1999), which was shown to return reliable turnover 
frequencies and be a promising candidate to account for 
the v p — LS anticorrelation, while relaxing considerably 
the constraints on the ambient conditions required by 
Bicknell et al.'s model. Yet following Begelman (1999), 
we ascribed FFA to the external layers of interstellar gas 
clouds that have been photoionized by the UV radia- 
tion from the active nucleus, and are actually ionization- 
bounded. These gas clouds might be associated with the 
NLR clouds. 

The particles of the lobes also produce inverse- 
Compton (IC) radiation via up-scattering of both the 



synchrotron radiation (synchrotron-self-Compton mech- 
anism; SSC) and the local, thermal photon fields. The 
energy density J7 ra d of the thermal fields was evaluated 
by taking into account the contributions by a putative ac- 
cretion disk, producing the bulk of its luminosity at UV 
frequencies, by a dusty torus, re-radiating the disk's UV 
photons in the IR domain, and by the stars of the host 
galaxy, mostly contributing the NIR-optical photons. We 
obtained C/yv ~ LS~ 2 and Um ~ LS~ 2 , whereas U opt is 
independent of LS for the considered LS < 1 kpc. 

The IC scattering of all the aforementioned radiation 
fields yields significant and complex high-energy emis- 
sion, from X-ray to 7-ray energies. In GPS quasars the 
putative direct X-ray emission of the accretion disk's hot 
corona and of the beamed relativistic jets may overcome 
the X-ray output of the lobes. In GPS galaxies these 
two contributions are instead expected, respectively, to 
be obscured by the torus and Doppler-hiddcn; the lobes 
are expected to be the dominant X-ray source. 

3. THE SOURCE SAMPLE 

We chose to apply our model to the sample of eleven 
GPS/CSO galaxies known as X-ray emitters up to 
2008. 12 The source list is given in Table Bl. 

The sources are all members of the catalogue of radio 
sources with flux density greater than 1 Jy at 5 GHz by 
Kiihr et al. (1981). The GPS nature of their radio spectra 
was observed at several epochs (O'Dea 1998, and refer- 
ences therein). Although the bona-ftde GPS nature of 
three of our sample's members (IERS B0108+388, PKS 
B0941-080, and IERS B1345+125) was recently ques- 
tioned on the basis of either flux-density variability or 
values of the radio spectral indices (Torniainen et al. 
2007), we decided to keep these sources in our sample be- 
cause the presence of the spectral peak at GHz frequen- 
cies was always confirmed. A CSO morphology charac- 
terizes the radio structures of all the sources, although 
one of them (PKS B1345+125) is an unusual example 
of CSO, with one of the two jets not being detected at 
frequencies above ~1 GHz. Our sample's members have 
radio power P 5GHz = 10 25 5 - 10 27 9 W Hz" 1 , and X-ray 
luminosities £2-iokcV — 6 x 10 41 — 5 x 10 44 erg s _1 ; they 
are hosted by galaxies located at redshifts z — 0.08—0.99, 
i.e. at luminosity distances between 0.4 and 6.6 Gpc. 
The linear (projected) sizes of the radio sources are in 
the range ^10 — 400 pc, thus providing "snapshots" of 
different stages of the source expansion. More details 
on classification and properties of individual GPS/CSOs 
and their host galaxies are given in Appendix A. The 
source physical quantities that are more relevant to our 
SED modeling are summarized in Tables Bl and B2. 

Fig. 1 shows the SEDs of the GPS/CSO galaxies of 
our sample. With the exception of some of the mid/far- 
infrared (MFIR) data, which will be discussed below, 
the SED data were derived from the literature (see Ap- 
pendix A for the complete reference list). In particu- 
lar, the near-infrared (NIR) and optical magnitudes were 
converted into fluxes by means of the absolute calibra- 
tions of Bessel (1979) for standard UBVK and Cousins' 

12 During the preparation of this manuscript, a paper by 
Tengstrand et al. (2009) appeared, reporting the detection of ad- 
ditional X-ray emitting GPS galaxies; the modeling of the broad- 
band spectra of these new sources is beyond the scope of the present 
work. 
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RI filters, Wamsteker (1981) for standard JHL bands, 
and Allen (1973) for standard RI filters; the dereddcning 
was performed by means of the extinction laws given by 
Rieke & Lebofsky (1985) and Cardelli et al. (1989), as- 
suming as B-band Galactic extinctions (Aq) the values 
by Schlegel et al. (1998) provided by the NASA/IPAC 
Extragalactic Database (NED). Source- intrinsic redden- 
ing effects were not considered in these bands, given the 
dominance of the host-galaxy contribution. The X-ray 
unabsorbed fluxes were derived from the absorbed fluxes 
by means of XSPEC v. 11.3.1, taking both the Galac- 
tic and the source-intrinsic hydrogen column densities 
(7VH,Gai and A^H.intr) into account. The SED luminosi- 
ties at the source rest-frame frequencies were then com- 
puted by assuming modern cosmology (see above). No 
if -correction was applied to the data. 

In our sample, literature data in the MFIR band were 
available for two sources only: IERS B1345+125 and 
IERS B1404+286. Given the importance of the MFIR 
portion of the SED for our modeling (see Section 4), we 
thus searched for and analysed archival Spitzer 13 IRAC 14 
(Fazio et al. 2004) and MIPS 15 (Rieke et al. 2004) obser- 
vations of our sample's members. 

IRAC data for IERS B1345+125 (prog. ID 32, January 
2004; Fazio & Wang 2004) and IERS B1404+286 (prog. 
ID 30443, July 2006; Rieke et al. 2006) were available in 
all of the four bands (3.6, 4.5, 5.8, and 8.0 fxm). MIPS 
observations at 24 and 70 fim were performed for IERS 
B1345+12 (prog. ID 30877, July 2007; Evans et al. 2006) 
and IERS B1404+286 (prog. ID 30443, July 2006; Rieke 
et al. 2006) an additional 24 ^m zodiacal light calibration 
image with IERS B0500+019 observed in the field (prog. 
ID 1882, April 2007) was also available. The analysis of 
the aforementioned IRAC and MIPS data is described in 
Appendix B; the flux densities are reported in Table B2. 

4. MODELING THE BROAD-BAND SPECTRA 

Our dynamical-radiative model, summarized in Sec- 
tion 2, is a powerful tool to study the evolution of a 
typical GPS-source synchrotron SED as a function of the 
time-dependent source linear size LS(t), given the kinetic 
power of the jets, Lj, and the energy spectrum of the 
hot-spot electrons injected into the lobes; furthermore, 
assuming typical luminosities for the putative torus' and 
accretion-disk's components, as well as for the host galax- 
ies, the model enables one to investigate the temporal 
evolution of the high-energy, comptonized SED compo- 
nent (Stawarz et al. 2008). 

By applying the model to sources with measured (pro- 
jected) linear sizes, we could here constrain their jet ki- 
netic powers, as well as the spectra of their hot-spot par- 
ticles; furthermore, we could test the viability of the FFA 
effect as the main responsible for the optically-thick part 
of the radio spectra; finally, the observational constraints 
on the luminosities of the host galaxies, and, for a few 
sources, of the torus emission, enabled us to evaluate the 
contribution of the comptonized radiation to the high- 
energy emission of the source. 

In Fig. 1, we show the modeling of the intrinsic broad- 
band SEDs of our sample's sources, and Table B3 reports 

13 http:/ /www. spitzer.caltech.edu/ 

14 Infrared Array Camera 

15 Multi-Band Imaging Photometer for Spitzer 



on the values of the best-fit parameters. More details on 
the modeling of the various SED components are given 
below. 

4.1. Synchrotron emission 

In our model, the computation of the theoretical radio 
spectrum is performed under the assumptions that the 
radio emission is dominated by the contribution of the 
lobes, and that each lobe contributes the same amount 
of radiation. Although the theoretical source morphol- 
ogy is necessarily much simpler than the actual source 
structure, the above approximations reasonably fit our 
sample of GPS/CSO galaxies. 

For the selected sources, we thus modelled the radio 
spectrum as synchrotron radiation produced by the elec- 
tron population of the lobes, N e ("f), which represents 
the evolution of the injected particle hot-spot popula- 
tion Q("f). Among the possible injection functions Q{^f) 
discussed in Section 2, we chose the broken power law 
Q{l) ~ l~ S i with s = si for 7 < 7; nt , and s = s 2 for 
7 > 7int, the break being fixed at 7; n t ~ 2 x 10 3 ; a similar 
spectrum was indeed inferred from the broad-band mod- 
eling of the hot-spot emission of the radio galaxy Cyg 
A, and interpreted as the result of a transition between 
two different acceleration mechanisms: the cyclotron res- 
onant absorption, and the diffusive shock acceleration 
(Stawarz et al. 2007, and references therein). The lower 
and higher Lorentz factors of the emitting particles were 
chosen as 7 m j n = 1 and 7 max = 10 5 . The emission of 
this particle population well reproduces the optically- 
thin part of the spectrum, i.e. the emission at frequencies 
higher than the turnover. 

FFA effects, as those generated by a spray of interstel- 
lar clouds engulfed by the expanding lobes and mingled 
with the synchrotron-emitting gas, are included in the 
model in order to compute the optically-thick portion 
of the spectrum, i.e. the emission at frequencies lower 
than the turnover. The value of the turnover frequency 
enters the model as an input parameter (f P) i n tr); the cal- 
culation of this value would indeed require the physical 
and geometrical modeling of the absorber, and thus the 
introduction of additional free parameters in the model 
(see Begelman 1999). 

The combination of the above synchrotron emission 
with the adopted absorption model returned a spectral 
shape of the form L v oc L v ^v~ 2 (with L„^ ua the lu- 
minosity unaffected by FFA), and satisfactorily fits the 
shape of most of the optically-thick spectra without the 
need of assuming source inhomogeneities. This feature of 
our scenario is not shared by the SSA scenario. SSA ef- 
fects would return an optically-thick spectrum with slope 
independent of the optically-thin slope, and equal to 5/2 
for a homogeneous source; to obtain a flatter spectrum 
an inhomogeneous source should be assumed. Neverthe- 
less, the issue of homogeneity can also play a role in our 
scenario. In fact, our absorption model significantly de- 
viates from the data in the case of IERS B2352+495, a 
source whose spectrum is much flatter than predicted at 
frequencies below the turnover; in this spectral region, 
the source structure is indeed more complex than we as- 
sumed (Araya et al. 2010, and references therein). 

The spectral fits described above enabled us to con- 
strain Lj, si, and s 2 (see Table B3). 
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Fig. 1. — Modeling of the source-frame radio spectra and SEDs of the GPS/CSO galaxies of our sample. Radio to X-ray data were 
derived from the literature (the references are listed, source by source, in Appendix A), excepted for the MFIR data of IERS B0500+019, 
IERS B1345+125, and IERS B1404+286, which were derived from the analysis of archival Spitzer data (Appendix B). The model curves, 
obtained with the parameters listed in Table B3, are displayed with different-style lines. Left panels: The solid-line curve shows the 
optically-thin synchrotron emission at frequencies higher than the turnover frequency, and the synchrotron free-free absorbed spectrum 
at frequencies below the turnover. Right panels: The black, solid lines represent the synchrotron emission; the black, dash-dotted line 
indicates the corresponding SSC emission; the grey, dashed line shows the thermal MFIR emission from the torus; the grey, dotted line 
shows the thermal star light; the black, dashed and dash-dot-dotted lines show the comptonized thermal emission from the torus and the 
disk, respectively; the black, dotted line indicates the comptonized star light. Note that the shown X-ray spectra are unabsorbed, whereas 
the radio spectra are displayed as absorbed. When an IC spectral component does not appear in the plot, its luminosity is below the scale 
minimum. 



4.2. Thermal radiation fields 

The thermal radiation fields that we considered as soft- 
photon sources for the IC emission are the host galaxy's 
stars, the circumnuclear torus, and the accretion disk. 

We are interested in evaluating the relevance of the 
contribution of these radiation fields in the IC process, 
rather than in carefully modeling the SED of these com- 
ponents. Therefore, for the sake of simplicity, we mod- 
elled the spectra of the above thermal photon fields as 
black-body spectra with the appropriate frequency peaks 
and bolometric luminosities. 

4.2.1. Host galaxy 

As mentioned in Section 1, the properties of GPS- 
source host galaxies are consistent with those of non- 



passively evolving giant ellipticals: they often display 
signs of recent star formation as well as morphological 
disturbances. Nevertheless, departures from our assump- 
tion of a constant-density medium for the inner kpc-scale 
core (see Section 2) are negligible for any GPS-source 
host galaxy. 

As shown by, e.g., Silva et al. (1998), the typical SED 
of an evolved giant elliptical displays a dominant optical- 
NIR hump peaking about 1-2 /im, with peak luminosity 
of a few ~10 44 erg s _1 , and a secondary MFIR hump 
peaking about ^100 yum, and a factor ^300 fainter. 
Starburst, interacting, and ULIRG galaxies can however 
show a MFIR hump with a peak luminosity up to a factor 
~50 higher than the optical-NIR peak. 

We could satisfactorily model all of our galaxy optical- 
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Fig. 1. — Continued 
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Fig. 1. — Continued 



NIR SEDs with a black-body spectrum peaked at A = 
1.5 fim (y = 2 x 10 14 Hz), and a bolometric luminosity 
in the range L opt = 3 x 10 44 - 10 45 erg s" 1 (see Table 
B3). MFIR data were available for three objects of our 



sample: in all of these cases, we did not associate the 
MFIR emission with the galaxy, but with circumnuclear 
dust (see Section 4.2.2). 
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4.2.2. Torus 

Besides the emission from the host-galaxy dust accom- 
panying possible star-formation activity, a putative cir- 
cumnuclear dusty torus re-radiating part of the UV ra- 
diation absorbed from the accretion-disk would also con- 
tribute to the total emission in the MFIR domain. 

Most of our sample's sources lack MFIR observations 
and/or detections. In the three objects detected at MFIR 
frequencies, the radio-to-MFIR spectra clearly show that 
the MFIR data lie well above the extrapolation of the 
synchrotron spectrum, ruling out the jet as the source 
of this radiation, and suggesting a thermal origin of the 
MFIR spectral component. 

These sources are not spatially resolved in this en- 
ergy window, preventing us to locate the site of the 
MFIR emission. Nevertheless, as we already mentioned 
in Stawarz et al. (2008), hints for the presence of a dusty 
torus around their AGN come from numerous observa- 
tions of samples of GPS and CSS sources. Heckman et 
al. (1994) showed, by means of IRAS data (12-100 /zm), 
that the MFIR luminosity of GPS and CSS sources is 
comparable to that of extended sources with similar ra- 
dio power and redshift, and is consistent with thermal 
emission from dust heated by the UV emission of the 
central AGN, rather than with radiation from a circum- 
nuclear starburst. Analysis of ISO data by Fanti et al. 
(2000) also proved that FIR (60-100 fim) luminosities of 
GPS and CSS radio galaxies are not significantly differ- 
ent from those of extended objects, and associated the 
observed luminosities (on average > 6 x 10 11 L Q ) to dust 
with temperatures from 25-30 K (warm component) up 
to 60-100 K (hot component), whose masses are ~ 5x 10 5 
M© and 2 x 10 s M Q respectively. Finally, combining ISO 
data at 2-200 /um with millimetric and sub-millimetric 
observations of 3CR radio galaxies and quasars, Haas 
et al. (2004) demonstrated that both classes show MIR- 
to-FIR luminosity ratios typical of powerful AGNs; the 
MFIR emission can be ascribed to dust with tempera- 
tures between ~30 and 100 K; quasars and galaxies pref- 
erentially display Lmir/^fir > 1 and < 1, respectively, 
according to the idea that the thermal dust emission is 
located in the torus, whose internal, hotter component 
cannot be observed by edge-on observers. The above 
findings were confirmed by Spitzer/MIPS (24 — 160 (im) 
observations of FRII radio galaxies and quasars: Shi et 
al. (2005) showed the dominance of thermal emission 
from AGN-heated dust in their MFIR spectra, and no 
difference between small- and large-scale sources; CSS 
sources with sizes of 2-400 kpc displayed MFIR lumi- 
nosities between 10 10 — 10 14 L , with no size-luminosity 
correlation. 

On the other hand, Pollctta et al. (2008) modelled, by 
means of clumpy-torus models, the IR SEDs of a sam- 
ple of MIR-luminous AGNs at z = 1.3 — 3 discovered 
by Spitzer (IRAC and MIPS, A = 3.6 - 160/zm): the 
MIR data could be reproduced by a single-temperature 
(^300 K) dust component; for sources detected at FIR 
frequencies as well, the additional, FIR spectral com- 
ponent could be modelled either with an AGN-heated 
colder dust component detached from the torus, or with 
an exceptionally powerful (L > 3.3 x 10 12 L Q ) starburst, 
thus leaving the issue of the origin of FIR emission open. 

Finally, Spitzer data of FRI radio galaxies (Leipski et 



al. 2009) recently proved that star formation is actually 
present in these sources, with a contribution to the MIR 
emission ranging from minor to dominant; for sources 
in which the nuclear MIR component could be identi- 
fied, this component is dominated by non-thermal emis- 
sion, although some FRIs do exist with clear indication of 
AGN-heated dust emission. The debate on MFIR emis- 
sion is thus clearly open. 

As discussed in Section 3, MFIR data were available 
for three of our sample's members only. For the sake 
of simplicity, we chose to model the MFIR emission of 
these GPS/CSO galaxies with single-temperature black- 
body spectra, which we associated with torus' dust. The 
spectral fit yielded black body temperatures T ~ 50 — 
125 K, and thus peak wavelengths Air = 23 — 60 ^m. 
For the remaining sources, with no MFIR data available, 
we assumed the torus' dust to radiate as a black body 
with temperatures T ~ 50 K; its spectrum thus peaks at 
wavelengths Air ~ 60 ^m. The broad-band spectral fit 
procedure yielded MFIR luminosities Lir = 8 x 10 43 — 
6 x 10 45 erg s" 1 (see Table B3). 

Even though our analysis of the Spitzer/IRAC ob- 
servations suggests the presence of a hotter, MIR- 
emitting dust component in two sample's sources (IERS 
B1345+125 and IERS B1404+286), these objects are 
brighter in the FIR than in the MIR domain. There- 
fore, our single-temperature black-body spectral model 
accounts for the dominant IR component. Furthermore, 
the exact shape of the IR spectrum does not play a cru- 
cial role in determining the resulting X-ray spectrum, 
which is more heavily affected by the peak frequency of 
the IR spectrum and by the electron energy distribution. 

4.2.3. Accretion disk 

Both the direct UV emission from the disk and the 
possible X-ray emission from a hot corona (e.g. Koratkar 
& Blaes 1999; Cao 2009) are expected to be largely ob- 
scured by the dusty torus, in edge-on AGNs. Emission 
lines are thus used to investigate the accretion properties 
of obscured sources. 

According to the findings by Baum & Heckman (1989) 
and Rawlings & Saunders (1991) for extended sources, 
the [Oni]A5700 luminosity was claimed to be positively 
correlated with the radio luminosity also for a sample of 
GPS and CSS sources (Labiano 2008): this would sup- 
port the traditional scenario according to which the AGN 
powers both the ionized gas and the radio emission. 

Additional sources of ionization might however be 
plausible. Early studies by Morganti et al. (1997) showed 
that no significant differences in line luminosities and ra- 
tios were present between compact (CSS) and extended 
radio sources with comparable radio power and redshift, 
and only tentative evidence for a lower Oni luminosity in 
CSS sources was found. More recently, however, Labiano 
(2008) claimed the existence of a possible positive corre- 
lation between Oin luminosity and the size of the radio 
source, suggesting that the jet expansion contributes to 
enhance the line emission. 

Kawakatu et al. (2009b) recently compared the ob- 
served properties of some narrow emission lines of a sam- 
ple of young radio-loud galaxies (average size of ^3 kpc) 
with those of a sample of radio-quiet Seyfert 2 galaxies, 
and found that young radio-loud galaxies have systemat- 
ically larger low- vs. high-ionization emission line ratios 
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([Oi]A6300/[Oiii]A5700). They concluded that power- 
ful, young AGNs favour accretion disks without a strong 
big blue bump (hereafter BBB), whereas a strong BBB 
would be present in Seyfert 2s. The BBB is actually 
missing in geometrically thick and optically thin accre- 
tion disks. However, such disks would be characterized 
by a radiatively inefficient accretion flow (RIAFs; e.g., 
Narayan & Yi 1995), unable to generate the bolometric 
luminosities necessary to produce the Olll emission ob- 
served in young radio galaxies (Kawakatu ct al. 2009b), 
and in conflict with the high accretion rates found by Wu 
(2009a) for a sample of young radio sources. 

As acknowledged by Kawakatu et al. (2009b) , an alter- 
native explanation for the observed line ratios might be 
a scenario in which young radio galaxies do have a disk 
SED with BBB, but they are characterized by system- 
atically different combinations of the hydrogen column 
density Ah and the ionization parameter U (generally, 
higher Ah and lower U). Alternatively, following our 
model, the hard spectral component needed to explain 
the Oi/Oiii line ratio might come from the UV photons 
produced through IC in the expanding lobes (see Sect. 
4.3). 

The above evidence highlights that the actual SED of 
the accretion disk in GPS sources is still not well con- 
strained, and might be different from that of both stan- 
dard disks and RIAFs. 

In our broad-band SED modeling, we assumed that the 
disks of GPS sources do produce the BBB (Koratkar & 
Blaes 1999) characteristic of Seyfert 2s and quasars: we 
fixed the emission peak at ~10 eV (i^jv = 2.45 x 10 15 
Hz; T ~ 1.16 x 10 5 K ) and found luminosities iuv = 
10 45 _ 1Q 46 ergs -i ( see Ta bie B3). This UV disk emis- 
sion would be partly reprocessed by the dusty torus and 
re-radiated in the MFIR band, partly absorbed by the 
NLR gas clouds, and partly upscattered off the lobe 
high-energy particles to 7-ray energies. In the modeling 
procedure, the disk luminosities were constrained by the 
MFIR luminosities (see Table B3); the MFIR luminosity 
values, in turn, were constrained either by MFIR obser- 
vations (when available) or by the X-ray spectral fitting 
(see Section 4.3). As we showed in Stawarz et al. (2008), 
possible detections of (or upper limits on) GeV 7-ray 
emission of GPS/CSO galaxies by the Large Area Tele- 
scope on board the Fermi Gamma-Ray Space Telescope 
(hereafter Fermi/LAT; Atwood et al. 2009) would con- 
strain the disk luminosity in the UV band, thus enabling 
a test of the BBB assumption for the disk spectrum (see 
Section 4.3). 

4.3. Inverse-Compton emission 

The high-energy spectral components were computed 
through inverse-comptonization of the aforementioned 
synchrotron and thermal radiation fields off the lobe 
electron population responsible for the synchrotron ra- 
dio spectrum. As discussed in Stawarz et al. (2008), we 
assumed a magnetic field only a factor of a few below 
the equipartition value (jib = 0.3, r] e = 3); this assump- 
tion is roughly consistent with the finding by Orienti & 
Dallacasa (2008b) that equipartition holds in the lobes 
of GPS sources. 

As Fig. 1 shows, in most of our sources the modelled X- 
ray emission is dominated by the inverse-comptonization 
of the MFIR photon field (IC(IR)). The contributions of 



the inverse-comptonized UV radiation (IC(UV)) and of 
the SSC emission reaches comparable magnitude in a few 
objects; the starlight does not play, as IC seed photon, a 
relevant role in any of our sources (see, however, Stawarz 
et al. 2008, for the complete theoretical case study). Our 
model well reproduces the observed X-ray spectra for 
reasonable assumptions on the thermal emission of the 
circumnuclear dust and of the accretion disk. 

The model also predicts significant 7-ray emission up 
to the GeV-TeV domain, not observed so far in any GPS 
galaxy. Even though, in some of our sample's sources, 
the high-energy tail of the IC(IR) spectral component 
reproducing the observed X-ray spectrum extends to the 
soft 7-ray energy band, the modelled 7-ray emission can 
be mostly identified as the result of the IC scattering 
of the UV photons produced by the accretion disk with 
characteristic BBB and appropriate luminosity, which we 
assumed to operate in the source. 

A comparison of our modelled SEDs with the 
Fermi/LAT sensitivity curves (Atwood et al. 2009) shows 
that none of our GPS/CSO galaxies could be detected 
by Fermi/LAT in less than one year: this is consis- 
tent with the current lack of any of these sources in 
the Fermi/LAT First Source Catalog (Abdo et al. 2010). 
The modeling also predicts that five out of eleven sam- 
ple's members (i.e., IERS 0026+346, IERS B0108+388, 
IERS B0500+019, IERS B1345+125, IERS B1404+286) 
could be detected by Fermi/LAT at 2a level after an in- 
tegration time of about one year. This prediction should 
however be taken with caution: as mentioned in Section 
4.2.3, the 7-ray fluxes depend upon the assumed accre- 
tion disk luminosities L\jy, as well as upon the high- 
energy tail of the lobe electron distribution. A non- 
detection by Fermi/LAT might thus imply either UV 
accretion disk luminosities lower than we assumed, or 
a different electron distribution in the lobes. 

5. FURTHER SUPPORT TO THE X-RAY LOBE SCENARIO: 
THE N H -N m CORRELATION 

5.1. Context 

As we recently pointed out (Ostorero et al. 2009), the 
prediction of the X-ray-emitting lobes that characterizes 
our model may be supported by further observational 
evidence. 

The X-ray emission of GPS galaxies was traditionally 
mostly interpreted as thermal radiation from the accre- 
tion disk, absorbed by a gas component associated with 
the AGN and characterized by an equivalent hydrogen 
column density Ah (O'Dea et al. 2000; Guainazzi et al. 
2004, 2006; Vink et al. 2006; Siemiginowska et al. 2008), 
rather than as non-thermal emission from the jet or the 
lobe. Non-thermal emission from the jet/lobe, although 
not always completely ruled out, was never considered to 
be significant. The main argument brought in support of 
a disk-dominated X-ray emission scenario is the apparent 
discrepancy between the equivalent total-hydrogen col- 
umn density Ah derived from the X-ray spectral analy- 
sis and the neutral hydrogen column density Ahi derived 
from the radio measurements of the 21-cm absorption 
line, when the latter is available. The authors concluded 
that, because Ah always exceeds Ahi of 1-2 orders of 
magnitudes, the X-rays must be produced in a source 
region that is more obscured than the region where the 
bulk of the radio emission comes from, and thus located 
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within the radio lobes; if this was not the case, an unrea- 
sonably high fraction of ionized hydrogen (Hll) should be 
assumed to account for the above difference (Guainazzi 
et al. 2006; Vink et al. 2006). Such a scenario would also 
be consistent with the observed anticorrelation between 
Ahi and linear size found by Pihlstrom et al. (2003), be- 
ing the fraction of ionized gas likely low in a young radio 
source with still expanding Stromgren sphere (Vink et 
al. 2006). 

The discrepancies between the Ah and Ahi values 
mentioned above should in fact be regarded with caution. 
The ./Vhi estimate is derived, from the measurements of 
the hydrogen 21-cm absorption lines, as a function of the 
ratio between the gas spin temperature T s and its cov- 
ering factor Cf, representing the fraction of the source 
covered by the Hi screen (e.g. Gupta et al. 2006). The 
common assumption T s /cf = 100 K refers to the case of 
complete coverage (c/ = 1) of the emitting source by a 
standard cold (Tk ~ 100 K) ISM cloud in thermal equi- 
librium, and thus with spin temperature T s equal to the 
kinetic temperature Tk (T s = Tk). However, this assump- 
tion returns a value of Ahi that represents a lower limit 
to the actual column density of both the total and the 
neutral hydrogen (e.g. Pihlstrom et al. 2003; Vermeulen 
et al. 2003; Gupta et al. 2006). In fact, when T k ~ 100 
K, the expected high abundance of the molecular hydro- 
gen, i?2 (Maloney et al. 1996), might account for the 
difference between Ah and Ahi . On the other hand, 
in the AGN environment, illumination by X-ray radia- 
tion might easily raise Tk up to 10 3 — 10 4 K (Maloney 
et al. 1994; Conway & Blanco 1995), making T s raise 
accordingly (Davies & Cummings 1975; Liszt 2001); a 
source covering factor smaller than unity would also in- 
crease the T s /cf ratio; both the above effects would po- 
tentially lead to Ahi values fully consistent with the Ah 
estimates. Finally, temperatures as high as several 10 3 K 
would likely imply the presence of a non-negligible frac- 
tion of ionized hydrogen (Hll; Maloney et al. 1996; Vink 
et al. 2006), also contributing to relax possible residual 
column-density discrepancies. 

Evidence is mounting that, in GPS and CSS sources, 
the Hi absorption lines are not generated by a screen 
covering the source uniformly: instead, they originate in 
clouds of neutral hydrogen connected with (at least one 
of) the jet/lobe radio structures, and possibly interacting 
with it (Morganti et al. 2004b; Labiano et al. 2006; Ver- 
meulen et al. 2006; Araya et al. 2010). The association of 
the bulk of the Hi absorption with the optical emission 
line currently supports the identification of the absorbers 
with the atomic cores of the NLR clouds, although the 
presence of Hi elsewhere is not ruled out (Labiano et al. 
2006; Vermeulen et al. 2006). 

In our GPS-source model, the NLR clouds are grad- 
ually engulfed by the expanding radio lobes. Moreover, 
the cloud density is assumed to decrease with the dis- 
tance r from the AGN as r~ n (with 1 < n < 2), according 
to the observations of the NLR in Seyfert galaxies (Krac- 
mer & Crenshaw 2000; Kraemer et al. 2000; Bradley et al. 
2004) , and consistently with the anticorrelation between 
Ahi and linear size found by Pihlstrom et al. (2003). 
Because in our model the contribution of the lobes to 
the X-ray output of the source is significant, becoming 
stronger and stronger with increasing source compact- 
ness, the NLR clouds might thus play an important role 
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Fig. 2. — X-ray column densities (TVjj ) as a function of ra- 
dio column densities (7Vjji ) for the seven GPS/CSO galaxies of 
our sample {sub-sample D5+U2; bullets), IERS B0108+388, IERS 
B0500+019, PKS B0941-080, IERS B1031+567, IERS B1345+125, 
IVS B 1358+624, and IERS B2352+495, and the three additional 
GPS galaxies reported by Tengstrand et al. (2009) (sub-sample 
ADU; diamonds): 4C +32.44, PKS 0428+20, and 4C +14.41. 
Black symbols: Nm was computed by assuming T B = 100 K and 
Cf = 1; arrows represents upper limits. Grey symbols: as an ex- 
ample, the same sources with T B /c f = 5.5 x 10 3 K. Dash-dotted 
line: linear fit to sub-sample D5+AD of Nn/Nm detections (with 
T s = 100 K); dotted line: linear fit to sub-sample D5+ U2+AD U in- 
cluding both detections and upper limits. Data are from: Guainazzi 
et al. (2006); Mirabel (1989); O'Dea et al. (2000); Pihlstrom ct al. 
(2003); Sicmiginowska et al. (2008); Sicmiginowska (priv. comm.); 
Tengstrand et al. (2009); Vermeulen et al. (2003); Vink et al. 
(2006). 



in the absorption of the lobe X-ray radiation. 

5.2. Ah -Ahi Correlation 

As discussed above, besides the modeling of the broad- 
band SEDs, a way of discriminating among different sce- 
narios, and unveil the actual X-ray production site, is 
to compare the properties of the X-ray and radio ab- 
sorbers, i.e. compare the column densities derived from 
the analysis of the X-ray data, Ah , and those obtained 
from radio measurements at the rcdshifted 21-cm wave- 
length, Ahi • Such a comparison can be performed on 
a source-by-source basis: an ad hoc increase of cither 
the T s parameter or the T s /cf ratio can always remove 
possible Ah and Ahi discrepancies. Alternatively, one 
can compare the column densities of a sample of sources: 
the existence of a positive, significant Ah-Ahi correla- 
tion would suggest that the X-ray and radio absorbers 
coincide, thus supporting the co-spatiality of the X-ray 
and radio source. 

We investigated the existence of a connection between 
Ah and Ahi in our sample. For a positive correlation, 
we searched the source sub-sample for which both Ah 
and ./Vhi estimates (either detections or upper limits) are 
available (see Table B4). The results of the correlation 
analysis are reported in Tables B5, B6, and Fig. 2, and 
are discussed below. 

Estimates of Ah were available for ten out of eleven 
sample's members, and an upper limit was available for 
the remaining source; Ahi was instead measured in 6 
sources only, and Ahi upper limits were derived for 2 
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additional objects. Estimates of both Ah and Ahi were 
available for a sub-sample of five sources; taking the up- 
per and lower limits to Ah and Ahi into account, the 
sub-sample extended to eight sources. 

Note that, due to the poor photon statistics, the up- 
per limit to Ah for source PKS B0941-080 was derived 
(A. Siemiginowska, priv. comm.) from a re-analysis of 
the Chandra data (Siemiginowska et al. 2008) in two dif- 
ferent ways: (i) by fixing the absorption equal to the 
Galactic value, and letting the photon index L free to 
vary; (ii) by letting the intrinsic absorption free to vary, 
and fixing the photon index to a value varying in the 
range L=1.7— 1.9, the upper-limit of this range being re- 
ported as the mean L of both radio galaxies (Brinkmann 
et al. 1995) and GPS sources (Brinkmann et al. 1997). 
In our analysis, we adopted the more conservative es- 
timate (i). For source IERS B1345+125, estimates of 
Ah were derived from two observations carried out by 
different instruments (ASCA and Chandra) in different 
epochs (1996 and 2000, respectively). Given the fact that 
we cannot rule out long-term column-density variations, 
we included both Ah estimates in our data set; however, 
including the average Ah does not change our results 
significantly. Finally, the physics of IERS B1404+286 
appears to be so complex that a variety of models can 
satisfactorily fit the X-ray data (Guainazzi et al. 2004). 
As a consequence, the inferred values of Ah can vary by 
a factor ^100 (see Table B4). This forced us to remove 
this source from our correlation analysis. 

Our final correlation sample thus consisted of seven 
sources. Hereafter, the sub-sample of objects for which 
column density estimates did exist will be referred to as 
sub-sample D5, wherease sub-sample U2 will indicate the 
sources for which upper limits to Ah and/or Ahi were 
estimated; our full correlation sample will be indicated as 
sub-sample D5+U2. The column-density detections and 
upper limits that entered the above sub-samples were 
indicated with boldface characters in Table B4. 

A Pearson correlation analysis applied to sub-sample 
D5 revealed a significant (probability of the null- 
hypothesis of no correlation being true: Prob. ~ 1.7 x 
10~ 5 ) Ah- Ahi positive correlation. The statistical signif- 
icance of this correlation, however, is admittedly driven 
by source IERS B0108+388, characterized by the high- 
est Ah-Ahi values. Indeed, in the same sub-sample, the 
evaluation of the correlation by means of more robust, 
non-parametric (or rank) methods (i.e. Spearman's and 
Kendall's correlation coefficients; Press et al. 1992), de- 
creased the significance of the correlation. 

In the full correlation sample, i.e., sub-sample D5+U2, 
the correlation was investigated by means of survival 
analysis techniques. In particular, we made use of the 
software package ASURV Rev. 1.2 (La Valley et al. 1992), 
which implements the methods for bivariate problems 
presented in Isobe et al. (1986). The sub-sample D5+U2 
was shown to display a positive correlation, with signifi- 
cance level about 8%. The results of our correlation tests 
are reported in the first raw of Tables B5 and B6. 

It is rather promising that the significance of the corre- 
lation substantially improved when the sample size was 
increased to ten sources by including the three new X-ray 
GPS galaxies, among those reported by Tengstrand et al. 
(2009), for which both Ah and Ahi estimates were avail- 



able. 16 In this case, the probability of the null-hypothesis 
of no correlation being true decreased to ^2%. The re- 
sults on the full composite sample, including the addi- 
tional sub-samples AD (additional detection) and ADU 
(additional detection and upper limits), are given in the 
second raw of Tables B5 and B6. We wish to stress that 
dropping source IERS B0108+388 (which drove the cor- 
relation in sub-samples D5 and D5+U2) from the full 
sample D5+U2+ADU, returned a correlation with still 
good significance (5%: Spearman; 6%: Kendall). 

We note that our tentative correlation, that we an- 
ticipated for a smaller source sample in Ostorero et al. 
(2009), is in agreement with the tentative anticorrela- 
tion between Ah and linear size reported by Tengstrand 
et al. (2009), given the known GPS-source Ahi -LS anti- 
correlation (Pihlstrom et al. 2003). We also note that a 
comparison of Ahi and Ah in a sample of spiral-hosted 
Seyfert galaxies did not reveal any correlation (Gallimorc 
et al. 1999). 

As far as the law describing the relationship between 
Ah and Ahi is concerned, according to Pearson's test 
we could fit a linear relation to the log(AH)-log(AHi ) 
sub-samples D5 and D5+AD; however, this relation is 
not a good description of the data (xj?ed = 9-63 and 9.80, 
respectively). As for the censored sub-samples, D5+U2 
and D5+U2+ADU, we applied the ASURV Schmitt's lin- 
ear regression to estimate the best-fit straight-line pa- 
rameters; however, this procedure did not enable us to 
evaluate the goodness of the fit. The results of our linear 
regression analysis are reported in the last two columns 
of Tables B5 and B6, and are displayed in Fig. 2. 

6. DISCUSSION 

Our systematic SED analysis provided us with con- 
straints on some interesting physical parameters of our 
sample's members, as the spectrum of the electrons in- 
jected from the hot spots into the lobes, the luminosities 
of the accretion disk and the torus, the accretion effi- 
ciency, and the jet kinetic power, as discussed below. 

The broken power-law energy distribution assumed for 
the electron population injected into the lobes is char- 
acterized by a break energy 7; nt ~ m p /m c ; furthermore, 
the lower-energy segment (7 < 7; nt ) of the spectra have 
indices 0.7 < si < 2, and are thus mostly flatter than the 
canonical spectra generated by diffusive (1st order Fermi) 
shock acceleration. These findings are in good agreement 
with the results inferred from the modeling of the broad- 
band spectra of both the hot spots of some extended 
radio galaxies (Stawarz et al. 2007; Godfrey et al. 2009) 
and luminous blazars (Sikora et al. 2009). Therefore, as 
the aforementioned cases, our sources seem to fit into a 
scenario in which two different acceleration processes are 
at work in mildly-relativistic shocks: a pre-acceleration 
process responsible for the lower-energy spectra, as e.g. 
the cyclotron resonant absorption (Hoshino et al. 1992; 
Amato et al. 2006) , and an acceleration process acting at 
higher energies, as e.g. the Fermi-type mechanism. The 
transition between the two acceleration regimes, marked 
by the energy break 7i nt , would reflect the dominant role 
of the protons in the jet's dynamics. 

The torus luminosities yielded by our modeling 

16 Detection: 4C +32.44; upper limits: PKS 0428+20 and 
4C +14.41 . 
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Fig. 3. — The Eddington ratios of the GPS/CSO galaxies of our 
sample, inferred from the SED model parameters, arc displayed 
as a function of their BH mass (data are from Table B3), and 
compared to the theoretical curves of the accretion disk instability 
model by Czerny et al. (2009). Source names are indicated in a 
short format in the plot. Filled circles indicate the sources for 
which m/m Edd was computed with the estimate of the BH mass; 
open circles show the sources (names are in parenthesis) for which 
the BH mass was assumed as the average BH mass of the GPS- 
source sample by Wu (2009b). Dotted and dashed-dotted lines 
represent the theoretical curves of constant duration (10 3 and 10 4 
years, respectively) for accretion disk outbursts triggered by the 
radiation pressure instability, assuming a disk viscosity a = 0.02 
(Czerny et al. 2009, their Fig. 4). 



(~10 44 — 10 45 erg s _1 ) are as large as the ones observed 
in powerful quasars and FRII radio galaxies (e.g., Shi 
et al. 2005); moreover, the accretion disk luminosities 
(10 45 — 10 46 erg s _1 ) are as large as the ones estimated 
for powerful fiat-spectrum radio quasars (FSRQs) (e.g., 
Sambruna et al. 2006; Ghisellini et al. 2009). However, 
the jet kinetic powers that we found (~ 2 x 10 44 — 4 x 10 45 
erg s _1 ) are smaller than the values claimed for FSRQs 
(10 46 -10 48 erg s" 1 ; e.g., Sambruna et al. 2006; Xu et al. 
2009; Ghisellini et al. 2009); they rather seem to be com- 
parable to the jet powers estimated for steep-spectrum 
radio quasars (SSRQs) with similar BH mass, and to be 
higher than those of BL Lacertae objects and FRI radio 
galaxies (e.g., Xu et al. 2009). 

As noted by Czerny et al. (2009), assuming a standard 
10% radiative efficiency for both the accretion disk and 
the radio jet, the jet kinetic power sets a lower limit to the 
accretion luminosity: i acc > Lj/10. The results of our 
modeling yielded accretion-disk luminosities Ltjv and jet 
kinetic powers Lj such that Ltjv > AjA^ (see Table B3, 
col. 10): this confirms the radiative efficiencies typically 
assumed, and characterizes our sources with jet /disk lu- 
minosity ratios of 0.01 — 0.1, a factor ^100 lower than the 
ratios typically found in powerful blazar sources (^1-10; 
e.g., Celotti & Ghisellini 2008). 

The best-fit values of the accretion-disk luminosities 
enabled us to make a reasonable guess on the source ac- 
cretion rates. BH mass estimates were available for eight 
of our sample's members; by using the estimated BH 
mass values for these sources, and the average BH mass 
of 1.698 x 10 8 Mq derived by Wu (2009b) for a sam- 
ple of GPS sources in the remaining three sources, we 
inferred for our GPS/CSO galaxies Eddington ratios of 



rh/riiEdd = Luv/LEdd = 0.026 — 0.506 (see Table B3, 
col. 11). 

These ratios correspond to accretion rates that are 
comparable to (for source IERS B1404+286) or much 
higher than (for all the other sources) the critical accre- 
tion rate rii* ~ 0.025 rriEdd proposed by Czerny et al. 
(2009) as the threshold above which the radiation pres- 
sure instability operates in accretion disks, giving rise to 
intermittent source activity. Therefore, the results of our 
modeling are consistent with all of our sample's members 
being intermittent radio sources, with IERS B1404+286 
being a border-line case. 

In this context, according to the prescription by Cz- 
erny et al. (2009), for values of the disk viscosity a in 
agreement with the observational constraints (a ~ 0.02; 
e.g., Starling et al. 2004), the sources characterized by 
m O.lrliEdd are potentially able to undergo outburst 
phases lasting ~10 4 years, and thus grow to super- 
galactic scales after a few of these outbursts, provided 
that they have the appropriate BH mass. Higher vis- 
cosity values would imply higher accretion rate thresh- 
olds. Except for IERS B1404+286, our sample's mem- 
bers are characterized by accretion rates higher than this 
limit; however, their BH masses imply outburst dura- 
tions shorter than 10 4 years. This can be seen in Fig. 
3, where our sample's members are plotted in the Mbh 
vs. rh/mEdd diagram, and are compared to the lines of 
constant outburst duration predicted by the radiation 
pressure disk instability model by Czerny et al. (2009): 
all the sources lie below the line representing the out- 
burst duration T = 10 4 years, and are thus expected to 
undergo shorter-lived outbursts which will prevent them 
to reach super-galactic scales. Interestingly, most of our 
sources follow the 10 3 year outburst model curve. Note 
that for three of our sample's members (open symbols in 
Fig. 3) the BH mass was assumed (see above); the avail- 
ability of a BH mass estimate might change the location 
of these objects in the diagram. 

7. CONCLUSIONS 

For the first time, we systematically compared the 
broad-band spectra of a sample of GPS/CSO galax- 
ies with a model. Specifically, we showed that the 
dynamical-radiative model we proposed in Stawarz et al. 
(2008) can reproduce the observed emission, from the 
radio to the X-ray energy range, of eleven GPS/CSO 
galaxies characterized by different linear sizes, and thus 
sampling different stages of the source expansion. 

The radio spectra were modeled as synchrotron emis- 
sion by a lobe electron population that represents the 
evolution of an injected hot-spot electron population suf- 
fering from adiabatic and radiative energy losses. At fre- 
quencies below the turnover, the shape of most of the 
spectra could be satisfactorily accounted for by an ab- 
sorption scenario dominated by the FFA mechanism. 

The X-ray emission could be interpreted as non- 
thermal radiation produced through IC scattering of the 
local thermal radiation fields off the lobe electron popu- 
lation. Among the possible seed photon fields, the dom- 
inant role in the generation of the X-ray radiation via 
IC scattering is played by the MFIR radiation, which 
we associated to the putative circumnuclear dusty torus; 
dust associated with circumnuclear star-forming regions 
might also contribute to this emission, but this would not 
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affect our results significantly. The UV radiation emit- 
ted by the accretion disk gives a contribution comparable 
to that of the MFIR photons in a few sources only, but 
provides the bulk of the seed photons responsible for the 
7-ray emission. Our model proved to be a viable alter- 
native to the thermal, accretion-disk dominated scenario 
for the interpretation of the X-ray emission of GPS/CSO 
galaxies. 

Finally, from a more observational perspective, our ra- 
diative model finds further support in the comparison 
of the measurements of the radio and X-ray hydrogen 
column densities of a sub-sample of our source ensem- 
ble: the data suggest a positive correlation, which, if 
confirmed, would point towards the co-spatiality of the 
radio and X-ray emission regions. New radio measure- 
ments, necessary to improve the statistics of our correla- 
tion sample, as well as the quality of the available data, 
arc currently being performed. 
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APPENDIX 

NOTES ON INDIVIDUAL SOURCES 
IERS B0026+346 (COINS J0111+3905; OB +343) 

The object, quite powerful in the radio band (Psghz = 10 27 1 W Hz -1 ), was observed repeatedly at many radio 
frequencies since the 1980's. Because of its convex radio spectrum (e.g., Jauncey et al. 1970; Kiihr et al. 1981), it was 
identified as a GPS source, although the spectral turnover did not look dramatic. The GPS nature of the source was 
also confirmed, on the basis of flux variability and spectral shape, by Torniainen et al. (2007). Despite the detection of 
a ~4. 7c apparent speed, atypical for GPS sources, for the brightest source component (Kellermann ct al. 2004), recent 
flux-density monitoring measurements by Torniainen et al. (2007) detected only moderate variability and a spectrum 
steep enough to confirm the classification of the object as a GPS source on the basis of their criteria. The source is 
seen as an unresolved core by the VLA (Ulvestad et al. 1981); however, VLBA maps enabled Kellermann et al. (1998) 
to classify the object as a double source, by assuming that the core of one of the weaker, intermediate components 
was the core of a CSO, although they could not exclude the core as being actually located at one of the extremities. 
The overall source structure extends over ~35 mas (~211 pc), and is hosted by a galaxy located at z = 0.517 (Zensus 
et al. 2002), whose optical i?-band magnitude is tor, ~ 21 (Snellen et al. 1996, and references therein). 

In the X-ray energy range, the source was firstly successfully detected (0.54-3.9 keV) by the Einstein Observatory in 
1980 (Brcgman et al. 1985; Kollgaard et al. 1995, and references therein). The target was then re-observed by ROSAT 
(0.1-2.4 keV), during both the ROSAT All Sky Survey (RASS) and a pointed observation in 1992 (Kollgaard et al. 
1995, and references therein), and was detected in a possibly lower flux level. Conclusive spectral information could 
be extracted neither from the Einstein nor from the ROSAT data. Finally, the source was the target of a pointed 
observation by XMM-Newton (0.5-10 keV) in 2004 (Guainazzi et al. 2006): the source was detected with a 2 — 10 keV 
luminosity of 2.3 x 10 44 erg s _1 . the data spectral fitting revealed a spectrum with slope a — 0.43^0 19 and a relatively 
high column density (Nh = 1.0~!i ' A x 10 22 cm~ 2 ), whereas the timing analysis showed the existence of a flux variable 
of a factor ~3 on time-scales of a few ks. From the comparison of the Einstein, ROSAT, and XMM 1-keV flux densities 
a soft -X-ray flux decline of a factor 3 over ~10 years emerged: however, due to the quality of the X-ray data, it was 
not possible to ascribe those variations to a change of the column density of the intervening absorber (Guainazzi et 
al. 2006). Based on the detection of rapid X-ray variability in the XMM data set, Guainazzi et al. (2006) constrained, 
through the light-travel argument, the size of the X-ray emitting region to iSlO /xpe, concluding that the X rays are 
produced either at the base of the jet or in the accretion disk; the lack of excess soft X-ray emission is also brought by 
the authors in support of a scenario in which the hotspots are X-ray silent. 
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References for spectrum and SED of Fig. 1 

Jauncey et al. (1970); Kapahi (1981); Ulvestad et al. (1981); Peacock et al. (1981); Perley (1982); Rudnick & Jones 
(1982); Hutchings et al. (1994); Heckman et al. (1994); Taylor et al. (1994); Kollgaard et al. (1995); Kovalev et al. 
(1999); Terasranta et al. (2001); Zensus et al. (2002); Hirabayashi et al. (2000); Waldram et al. (2003); Terasranta et 
al. (2005); Beichman et al. (1981); Stickel et al. (1994); Snellen et al. (1996); Guainazzi et al. (2006); NED. 

IERS B0108+388 (COINS 30029+3^56; OC +314) 

This GHz-peaked radio source, very powerful in the radio band (Aghz = 10 27 3 W Hz" 1 ), has been classified as a 
CSO. It displays a nucleus with inverted spectrum, and twin jets with steeper spectra extending 3 mas (—20 pc) from 
the nucleus toward the NE and SW (Carilli et al. 1998, and references therein). High-resolution imaging of the source 
was first performed at 5 GHz as part of the PR VLBI survey (Pearson & Readhead 1988), showing the source to have 
a simple double structure. Conway et al. (1994) and Taylor et al. (1996), with more sensitive and higher resolution 
images at higher frequencies, revealed a steeply inverted core connected to both outer components by a faint chain of 
components. Overall, IERS B0108+388 has an "S" symmetry, a relatively common feature among CSOs (Readhead 
et al. 1996). Unlike most CSOs, IERS B0108+388 also displays a large-scale jet with steep spectrum, extending about 
25" (—170 kpc) to the east of the nucleus (Baum et al. 1990). Proper motions enable an estimate of the kinematic age 
of -400 yr (Owsianik et al. 1998; Polatidis & Conway 2003). The radio flux density is weakly polarized (0.30 ± 0.08 
at 4.8 GHz); although it was claimed not to vary significantly (Aller et al. 1992), a study by Torniaincn et al. (2007) 
revealed flux variations greater than a factor of —6 at some radio frequencies, and led the authors to not confirm the 
true GPS nature of the source. The radio source is associated with a narrow emission line galaxy at z=0. 66847 with 
R=22.0 mag (Stanghellini et al. 1993; Stickel et al. 1996a; Carilli et al. 1998), and is located (in projection) within 
0.5" (—3.4 kpc) of the galactic center. Optical images of the host reveal a very red, diffuse, and slightly asymmetric 
galaxy (perhaps a face-on spiral), although the faintness of the galaxy makes the classification difficult (Stanghellini et 
al. 1993). There is no evidence for a strong point-source contribution in the R-band image. The I-band image is more 
compact, leading Stanghellini et al. (1993) to propose an increased contribution from the active nucleus in the NIR, 
possibly indicating the existence of nuclear obscuration toward the AGN in IERS B0108+388. Variability in the near 
IR, with a very red observed color (aiR^ opt < —3) during an IR maximum, was revealed by Stickel et al. (1996a). 

By means of the M R - M B h relation, Wu (2009b) estimated a BH mass of 7.94 x 1O 7 M . 

In the X-ray band, the source was detected for the first time by XMM-Newton in 2004 (Vink et al. 2006), with 
a luminosity of 1.18 x 10 44 erg s _1 . The source counts were not sufficient to independently estimate, in the fitting 
procedure, spectral index a and column density Nn . By fixing a = 0.75, Nn= (57 ± 20) x 10 22 cm" 2 was obtained. 

The detection of a strong Hi 21-cm absorption line, with a width of 100 km s" 1 and optical depth r = 0.44, yields a 
neutral hydrogen column density Nm— 80.7 x 10 20 cm -2 , under the standard assumption of T s — 100 K and c/ = 1 
(Carilli et al. 1998). VLBI opacity maps recently allowed to prove that the mechanism responsible for the turnover of 
the radio spectrum is SSA in the core, whereas in the lobes FFA by ionized gas with a temperature T — 2000 K (Marr 
et al. 2001). 

References for spectrum and SED of Fig. 1 
Tinti et al. (2005); Dallacasa et al. (2000); Stickel et al. (1996a); Vink et al. (2006); NED. 

IERS B0500+019 (PKS 0500+019; OG +003) 

Observed in the radio band since the 1970's, and well known as a convex-spectrum radio source (e.g., Jauncey et al. 
1970), this very powerful (P 5G Hz = 10 273 W Hz" 1 ) object was classified as a GPS source by O'Dea et al. (1991). The 
GPS nature of the source was recently confirmed, on the basis of flux variability and spectral shape, by Torniainen et 
al. (2007). VLBI, VLB A and VSOP images show a radio structure displaying a symmetric, S-shaped morphology, as 
other GPS sources do, dominated by components similar to jets and/or micro- lobes, with the possible presence of a 
weak core, and with the northern part brighter than the southern part (Stanghellini et al. 1997, 2001; Fey & Chariot 
2000; Fomalont et al. 2000). The total extension of the radio sources is —15 mas (—96 pc). The total fractional 
polarization is below 0.1% at 5 GHz (Stanghellini et al. 1998), although polarized emission at the — 0.3% level is 
detected in the brightest component (Stanghellini et al. 2001). 

The optical host of the radio source was the subject of a long debate: originally identified with an unresolved QSO 
lying 2" north of a galaxy with m r = 21.5 and with extremely red optical spectrum (a^B = 7.1, with F v oc ^" Q ; 
Fugmann et al. 1988; Fugmann & Meisenheimer 1988), it was subsequently associated with the southern, asymmetric 
galaxy, with total magnitude = 20.7 and dominated by a strong point-like source in the NIR band (Stickel et al. 
1996a). Because this galaxy is part of an apparent group of more than 10 galaxies, the observed asymmetry might be 
generated by gravitational interactions with companion objects (Stickel et al. 1996b). The galaxy produces moderate- 
ionization narrow emission lines (Stickel et al. 1996b), and is located at z = 0.583 (de Vries et al. 1995). Based on the 
detection of a further, unidentified emission line, Stickel et al. (1996b) proposed that a background, reddened quasar 
is strongly aligned with the foreground galaxy. This scenario, however, was not confirmed by subsequent observations 
by de Vries et al. (1998a, b, 2000). Besides confirming the galaxy at z = 0.583 to be the host of the radio source, these 
authors revealed a degree of source nucleation much higher in the K band than in the J and H bands, contrary to the 
behavior of both the GPS sources of the same redshift and the FRIIs with the same degree of J-band nucleation; this 
led the authors to suggest that the very red color of the nucleus is produced by an absorption excess in the host galaxy. 
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Jackson et al. (2002) too identified the host as a narrow-line radio galaxy, providing a revised redshift of z = 0.584 
(Hook et al. 2003). 

Whereas the radio emission of the source is fairly stable (Stanghellini et al. 1997; Torniainen et al. 2007), moderate 
variability was detected in both the optical and NIR bands (Stickel et al. 1996a). 

The source is also a MFIR emitter: we detected it in the field of an archival Spitzer/MlPS calibration observation 
at 24 (Um (see Section 3, and Appendix B for details). 

In the X-ray regime, the Einstein Observatory observed the source in 1980, yielding however only a flux upper limit 
in the 0.33 - 4.64 keV energy range (Ledden & O'Dell 1983). 

The source was successfully detected by ROSA T during the RASS as well as in two pointed observations in 1993: 
although a unique spectral fit was not possible due to the paucity of low-energy counts, the spectrum was found 
to be consistent with either a significant excess absorption or a hard slope, even though thermal spectral models 
(however with poorly constrained temperature) could not be ruled out (Kollgaard et al. 1995, and refs. therein). The 
source spectrum was better determined through the XMM-Newton observations in 2004 (Guainazzi et al. 2006): with 
a 2 — 10 keV luminosity of 5 x 10 44 erg s _1 , the source spectrum was well fitted with a slope a = 0.62±Q; 2 g and 

A H = 5.0t|o x 1021 cm ~ 2 • 

An independent estimate of the column density of neutral hydrogen comes from the measurements of the redshifted 
21-cm absorption lines, yielding Ahi= 6.2 x 10 18 (T S //) cm -2 (Carilli et al. 1998). Based on the standard assumption 
of T s ~ 100 K and c/ = 1, yielding Ahi= 6.2 x 10 20 cm~ 2 , Guainazzi et al. (2006) conclude that a difference of a factor 
~10 between Ahi and Ah supports a source scenario in which the absorbed X-rays come from a region located well 
within the lobe radio hotspots, completely X-ray silent. Assuming a non-thermal radio-to-optical quasar spectrum 
for the source, Carilli et al. (1998) derived, from the difference between the optical data and the extrapolation of the 
radio-to-infrared spectrum into the optical band, a lower limit to the rest-frame visual extinction, Ay > 3, implying 
T s > 500 K, and thus Ahi> 3.1 x 10 21 // cm~ 2 . This limit would be largely consistent with the above-mentioned Ah 
value provided by the X-ray data analysis. 

References for spectrum and SED of Fig. 1 

de Vries et al. (1998b); Stickel et al. (1996a); de Vries et al. (1995); Fugmann & Meisenheimer (1988); Guainazzi et 
al. (2006); this work (Spitzer/MIPS data); NED. 

IERS B071 0+439 (COINS J0713+4349; OI +417) 

This source is a well known CSO (Taylor et al. 1996), with radio power P 5G Hz = 10 27 1 W Hz" 1 . The flux density 
is very weakly polarised (< 0.15% ± 0.11% at 5 GHz), and the observed variations of the flux are not statistically 
significant (Allcr et al. 1992; Torniainen et al. 2007). The known GPS nature of the source was recently confirmed, on 
the basis of flux variability and spectral properties, by Torniainen et al. (2007). 

The total angular size of the source is 24.1 mas (Owsianik et al. 1998), which corresponds to a projected linear 
size of ^145 pc. VLBI maps at several frequencies (1.6, 5, 10.7, and 15 GHz) showed the overall triple structure of 
the source. Despite having three components, this source was provisionally classified as a compact double, based on 
the fact that more than 80% of the emission came from two almost equally bright components (Pearson & Readhead 
1988). Conway et al. (1992) argued that the two outer components were hotspots and minilobes, whereas the centre 
of activity was associated with the middle component, based on its compactness, spectrum and weak flux density 
variability. However, more recent multi-frequency observations reveal a compact component with a strongly inverted 
spectrum at the southern end of the middle component, suggesting that the true centre of activity lies there (Taylor 
et al. 1996). 

Readhead et al. (1996) gives age estimates of 1200-1800 yr based on synchrotron ageing, and 1500-7500 yrs from 
energy supply arguments; more recently, the first detection of hotspot advance velocities in a CSO was measured in 
this source by Owsianik & Conway (1998); these measurements, together with those by Polatidis & Conway (2003), 
yielded an age smaller than 1000 yr. 

In the optical band, IERS B0710+439 has been identified with a galaxy (Peacock & Wall 1981) of magnitude 
r = 19.7 ± 0.2 (Peacock & Wall 1981). The galaxy displays a very irregular morphology, with large differences in the 
two bands, and has a very red color (r — i = 1.9; Stanghellini et al. 1993). There is evidence of strong interaction, and 
presence of a large amount of obscuring matter. The numerous objects around the source suggest that 0710+439 is the 
dominant galaxy of a cluster of galaxies (Stanghellini et al. 1993). The emission-line redshift of the galaxy is z=0.518 
(Lawrence et al. 1996); the optical spectrum shows absorption lines characteristic of an evolved stellar population, and 
an optical continuum shape typical of an elliptical galaxy without any evidence for a non-stellar component. 

By means of the M R - M B h relation, Wu (2009b) estimated a BH mass of 2.512 x 1O 8 M . 

At high energies, the source was detected as an X-ray emitter by XMM-Newton in 2004 (Vink et al. 2006). Its 
2 — 10 keV luminosity is 2.16 x 10 44 erg s _1 , and the spectral analysis yielded a slope a = 0.59 ± 0.06 for the energy 
spectrum, and a source-intrinsic column density Ah= (0.44 ± 0.8) x 10 22 cm~ 2 . 

References for spectrum and SED of Fig. 1 
Stanghellini et al. (1993); Snellen et al. (1996); O'Dea et al. (1996); de Vries et al. (1998a); Vink et al. (2006); NED. 
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PKS B 094 1-080 

The source has a radio power Psghz = 10 26 1 W Hz -1 , and an overall size of ~50 mas (~ 177 pc) (O'Dca 1998). 
Imaging with the VLA shows a slightly resolved secondary component 20" (~71 kpc) east of the main one (Stanghellini 
et al. 1998). The VLBI morphology of the main component is that of a compact double (Dallacasa et al. 1998) typical 
of many CSS/GPS radio galaxies. The source was classified as a GPS source (O'Dea 1998, and references therein), 
although recent investigation by Torniainen et al. (2007) did not confirm the GPS-source classification, indicating the 
object as a steep-spectrum source. 

Ground-based optical-NIR observations showed the host galaxy of the main radio source as an elliptical envelope, 
with magnitude r = 17.9 and color r - i = 1, including two cores (Stanghellini et al. 1993). HST-NICMOS (JK) 
imaging revealed that the system is actually composed of two galaxies, interacting with each other, and characterized 
by blue nuclei (de Vries et al. 2000). The GPS source is associated with the core of the northern, larger, and brighter 
galaxy (Stanghellini et al. 1993; de Vries et al. 2000). The prominent emission line spectrum allows one to derive, from 
seven identified spectral features, a redshift z = 0.228. 

By means of the Mbh — & relation, Wu (2009a) recently estimated a central BH mass of 5.01 x 1O 7 M0. 

In the X-ray regime, the source was pointed at by Chandra in 2002 (Guainazzi et al. 2006; Siemiginowska et al. 
2008): the target was detected, and subsequent, independent analyses first failed in determining a significant spectral fit 
(Guainazzi et al. 2006) (fixed a — 1.0; Ah undetermined), and then derived a steep-spectrum source (a — 1.62~^\^) 
not absorbed by any excess hydrogen column (Siemiginowska et al. 2008). The latter analysis also reported the 
detection of the eastern source component revealed by the VLA, although with very low source counts. The 2-10 keV 
X-ray luminosity of PKS B0941-080 is the lowest of our GPS-galaxy sample (9 x 10 41 erg s _1 ; Guainazzi et al. 2006). 

References for spectrum and SED of Fig. 1 
de Vries et al. (1998a); Stanghellini et al. (1993); Guainazzi et al. (2006); Siemiginowska et al. (2008); NED. 

IERS B1031+567 (COINS 31035+5628; OL +553) 

This source, with radio power Psghz = 10 26 9 W Hz -1 , is a known GPS source (O'Dea 1998, and references therein), 
and was also included in the sample of bona- fide GPS sources assembled by Torniainen et al. (2007). 

The source was found to have a compact double morphology in 5 GHz observations by Pearson & Readhead (Taylor et 
al. 1996). VLBA observations at 8.4 and 15.4 GHz showed the source to be dominated by two leading edge-brightened, 
steep-spectrum components (Taylor et al. 1996), proposed by the authors to be the working surfaces and lobes of two 
oppositely directed jets. Based on these data, the authors concluded that the source was likely to be a CSO, with the 
center of the activity located in the source center. The CSO morphology was later confirmed by VSOP observations 
(Fomalont et al. 2000). The overall source size is ~40 mas (^224 pc). 

Measurements of the hotspot separation velocity, based however only on two epochs, yield estimates of ^600-1800 
yr (Taylor et al. 2000; Polatidis & Conway 2003). 

In the optical band, the source is identified with an elliptical-like galaxy with magnitude r = 20.2, showing a possible 
asymmetric morphology at low brightness levels; the presence of weak resolved objects in the field suggest that this 
galaxy is the dominant member of a galaxy cluster. 

By means of the Mr - M B h relation, Wu (2009b) estimated a BH mass of 1.259 x 10 8 M©. 

In the X-ray band, the source was detected for the first time in 2004 by XMM-Newton (Vink et al. 2006), with a 
2-10 keV luminosity of 2.2 x 10 43 erg s" 1 . The source counts were not sufficient to independently estimate, in the 
fitting procedure, spectral index a and column density Ah . By fixing a = 0.75, Ah= (0.5 ± 0.18) x 10 22 cm~ 2 was 
obtained. 

References for spectrum and SED of Fig. 1 
Stanghellini et al. (1993),0'Dca et al. (1996); Snellen et al. (1996); de Vries et al. (1998a); Vink et al. (2006); NED. 

IERS B 13^5+ 125 (PKS 1345+12; 4C +12.50) 

This is a luminous, compact (~170 pc) radio source hosted by an elliptical-like galaxy at z = 0.12174, classified as 
an ultra- luminous infrared galaxy (ULIRG). 

The main host-galaxy body displays the same shape from NIR to near-UV wavelengths, is ~15" (^32 kpc) wide 
and rather asymmetric, and appears to be associated with a group of 15 fainter nebulous objects, many of which are 
located very close to its distorted halo (Heckman et al. 1986). Ground-based and HST observations clearly reveal two 
nuclei separated by ~1.8" (^3.8 kpc) embedded in a common, distorted envelope, including a strongly curved tidal 
tail (e.g., Gilmore & Shaw 1986; Heckman et al. 1986; Smith & Heckman 1989; Labiano et al. 2008). 

The consistency of the inner (10") NIR brightness profile with that of a merger remnant, and the little evidence of 
extinction, indicative of ISM stripping, in the very central region (Shaw et al. 1992), together with the detection of a 
rich young stellar population (YSP) (Zaurin et al. 2007), confirmed that this complex has a merger origin, as firstly 
proposed by Gilmore & Shaw (1986) to explain the double nucleus and the asymmetry of the large-scale structure. 

The south-eastern (SE) nucleus is believed to belong to an elliptical galaxy (Gilmore & Shaw 1986), and its NIR 
color is consistent with reddened starlight (Evans et al. 1999). The north-western (NW) nucleus displays a narrow- 
line Seyfert-2 optical spectrum (Gilmore & Shaw 1986), although the presence of a buried quasar is suggested by 
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the detection of NIR broad Paa lines (Veilleux et al. 1997, 1999), by the similarity of the extreme NIR color with 
that of quasars with a warm (500—1000 K) dust component (Evans et al. 1999; Surace & Sanders 1999), and by the 
highly-polarized UV emission (Hurt et al. 1999); stellar velocity dispersion derived from NIR VLT spectra indicate the 
presence of a 6.54 x 10 7 M Q BH (Dasyra et al. 2006). 

The source, detected by IRAS and ISO as a bright MFIR emitter (Golombek et al. 1988; Fanti et al. 2000), was 
recently observed in the MFIR domain by Spitzer, with both IRAC and MIPS (see Section 3, and Appendix B for 
details). 

Although the optical counterpart of the radio source was originally identified with the SE nucleus (Gilmore & Shaw 
1986), newer astrometry led to the association of the radio activity to the NW nucleus (Stanghellini et al. 1993; Evans 
et al. 1999; Axon et al. 2000; Batchcldor et al. 2007). The latter identification is adopted throughout this paper. 

The powerful (P 5G Hz = 10 25 94 W Hz" 1 ; O'Dea 1998) radio source is resolved neither by MERLIN (at 18 cm) nor 
by VLA (at 6 and 2 cm) observations by Spencer et al. (1989), which constrained the (projected) linear size to < 0.1" 
(< 213 pc). VLBI images at 5 GHz show a single jet disrupting at ~30 mas (^60 pc) south of the core, and then 
expanding in a diffuse ~40-mas-sized lobe reaching a distance of ~80 mas (~170 pc) from the core (Stanghellini et al. 
1997; Shaw et al. 1992). The puzzling lack of any counter-lobe at these frequencies (Stanghellini et al. 1997, 2001), 
together with the observation of highly-polarized components and superluminal motion (Lister et al. 2003) make the 
source an unusual example of CSO. A more symmetric structure is however detected by the VLBI at lower frequencies 
(~ 1266 MHz), where a fainter counter-jet/lobe extending toward the North for 2240 mas can be identified (Morganti 
et al. 2004b). The radio structure is misaligned by ^45deg with respect to the UV light distribution (Axon et al. 
2000). 

The source was shown to be non- variable at 318 and 430 MHz over a 14-year period of flux-density monitoring with 
the Arecibo telescope (Salgado et al. 1999). Although classified as a GPS source (O'Dea 1998, and references therein), 
a recent study by Torniainen et al. (2007) did not confirm the GPS nature of the object, classifying the source as a 
steep-spectrum source. 

In the X-ray domain, IERS B1345+125 was first revealed by ASCA (2-10 keV) in 1996 (O'Dea et al. 2000), after a 
previous non-detection by ROSAT (0.2-2 keV; O'Dea et al. 1996), and was the first GPS source to be reported as an X- 
ray emitter. The ASCA image does not resolve the two nuclei. However, a comparison of the X-ray power with the Ha 
emission-line luminosity makes O'Dea et al. (2000) conclude that the X rays come from the western nucleus, associated 
with the GPS source. Highly uncertain spectral index (a = 0.6+0;!) an d column density (Vh= 4. 2+2 ° x 10 22 cm~ 2 ) 
are derived from the modeling of these X-ray data. Higher-quality X-ray data were obtained in 2000 with Chandra 
(Sicmiginowska et al. 2008), which also revealed the presence of extended emission on ~10 arcsec (20 kpc) scale. 

The detection of CO(l— >0) emission implies a mass of molecular hydrogen M(H 2 ) ~ (3.3 - 6.5) x 1O 1O M (Mirabel 
et al. 1989; Evans et al. 1999), concentrated in the central ~2 kpc of the galaxy. 

Hi absorption was detected against the radio source by Mirabel (1989), showing a narrower (FWHM~200 kms -1 ) 
and a broader (FWHM~700 kms -1 ) component, and implying an average Vhi=6.2 x 10 18 (T S //) cm~ 2 . More recent 
measurements by Morganti et al. (2004a, b) revealed a much wider (~2000 kms -1 ) broad line, blucshifted with respect 
to the systemic velocity of the galaxy. Morganti et al. (2004b) also assessed, thanks to the spatial resolution of the 
VLBI, the actual location of the Hi absorber, showing that a single Hi cloud of (projected) size <L 22 x 65 pc, with mass 
M ~ (10 5 - 10 6 )M Q , column density V H i~ 10 22 (T s /100iT) cur 2 , and located at the edge of the northern jet/lobe 
(projected distance of 40-100 pc from the core), can account for the whole unresolved Hi absorption. This would prove 
the clumpy nature of the ISM in the nuclear region of this source, making the scenario of a strong jet-cloud interaction 
a viable interpretation of the observations of not only this galaxy, but also of other young radio sources (Morganti et 
al. 2004b). Holt et al. (2003) found complex optical emission line profiles at the position of the nucleus with line width 
of ^2000 kms -1 and blueshifted with respect to the narrow HI absorption component. 

Whether the GPS source is truly young is still a matter of debate. A model-dependent estimate by O'Dea et al. 
(2000) implies the source to be younger than 3 x 10 4 yr. However, Zaurfn et al. (2007) find no clear evidence of a 
time delay of the AGN activity with respect to the major merger-induced starburst (<6 Myr ago); this would also be 
in agreement with the apparent connection between the extended (SJ100 kpc; Stanghellini et al. 2005) and compact 
radio emission, indicating recursive AGN activity on a time-scale longer than previously thought. Lister et al. (2003) 
speculate the existence of a precessing BH to explain the S-shaped jet; if confirmed, this would point to an age of the 
radio source <10 5 yr. 

References for spectrum and SED of Fig. 1 

Pushkarev et al. (2005); Stanghellini et al. (1993); Snellen et al. (1996); O'Dea et al. (1991); Golombek et al. (1988); 
dc Vries et al. (1998a); Fanti et al. (2000); Scoville et al. (2000); O'Dea et al. (2000); Siemiginowska et al. (2008); this 
work (Spitzer/IRAC, and Spitzer/MIPS data); NED. 

IV S B1358+624 (COINS J1400+6210; 4C +62.22) 

This source has radio power Psghz = 10 26 9 W Hz -1 ; it is one of the 11 low- frequency variable sources mapped with 
global VLBI at 608 MHz by Padrielli et al. (1991), who showed the object to be a compact double separated by 46 
mas with a possible bridge of emission. Dallacasa et al. (1995) confirm that the two lobes are connected by a knotty, 
broad jet; Taylor et al. (1996) provide new evidence of the source structure: they resolve the two lobes, and detect a 
long, one-sided jet with a faint, compact component at the base; based on the compactness of the component, on its 
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inverted spectrum, and on its location, the authors confirmed the source to be a member of the CSO class, with the 
central component being the center of the activity. No evidence of hotspots on either sides of the central engine was 
found (Taylor et al. 1996). The overall source size is ~70 mas (^380 pc). 

The convex, stable, and GHz-peaked spectrum enabled the classification of the source as a GPS source (O'Dea 1998; 
Torniainen et al. 2007). 

An optical image of this source (Stanghellini et al. 1993) shows that the host is a galaxy of r = 19.8 mag with a blue 
color (r — i = 0.3) and a boxy morphology; there is evidence of obscuring material in a direction (NE-SW) roughly 
perpendicular to the VLBI axis, but on larger scale (Dallacasa et al. 1995). 

In the X-ray band, XMM-Newton detected the source in 2004 as an X-ray emitter with 2 — 10 keV luminosity 
of 1.67 x 10 44 erg s -1 (Vink et al. 2006). The fitting procedure yielded a spectral index a = 0.24 ± 0.17 and a 
total-hydrogen column density jV h = (3.0 ± 0.7) x 10 22 cm~ 2 . 

By means of the Mr — Mbh relation, Wu (2009b) estimated a BH mass of 1.585 x 10 8 M Q . 

An independent estimate of the amount of neutral hydrogen comes from the 21-cm observations by Vermeulen et 
al. (2003), who detected a column density Vhi= 1-88 x 10 20 cm~ 2 , under the usual assumption of T s = 100 K and 
uniformly covered source. 

References for spectrum and SED of Fig. 1 
de Vries et al. (1998a); Stanghellini et al. (1993); O'Dea et al. (1996); Snellen et al. (1996); Vink et al. (2006); NED. 

IERS B1404 +286 (MRK 0668; OQ +208) 

Observed in the radio band since the 1970's, this compact radio source is one of the closest bright (P5GHZ = 10 25 ' 4 
W Hz -1 ) GPS galaxies. Its radio spectrum is convex, stable, and GHz-peaked, and the classification of the object as 
a GPS source (O'Dea 1998) was confirmed also by Torniainen et al. (2007). 

The radio source displays a CSO morphology (Stanghellini et al. 2002): it has a weak core at 15 GHz, two-sided 
faint and short jets, and two mini-lobes located 10 mas apart, giving a projected size of less than 10 pc. The two 
micro hot-spots of this radio source seem to increase their distance at a rate of about 0.3c, comparable to what already 
found in other CSOs, yielding a kinematical age of a few centuries (Stanghellini et al. 2002). 

Based on VSOP observations, Kameno et al. (2000) reported a very inverted spectrum below the peak frequency of 
the weaker SW region of this source. They considered the extremely inverted spectrum inconsistent with SSA, and 
propose a process of FFA to explain the radio spectral shape. Xiang et al. (2002) confirmed the incompatibility of the 
spectrum with pure SSA result, on the basis of ground-based VLBI observations. 

The host galaxy of the radio source, located at z = 0.077, has mji — 14.6 and color r — i = 0.2; companions in 
the galactic envelope, and a tail of low brightness emission in N-S direction, suggest that the galaxy is dynamically 
disturbed (Stanghellini et al. 1993). The optical spectrum is that of a Seyfert 1; however, given the power of the 
associated radio source, it was classified as a broad line radio galaxy (e.g., Stanghellini et al. 1997, and references 
therein) . 

IRAS detected IERS B1404+286 as a bright MFIR emitter (Golombek et al. 1988); recently, the target was re- 
observed in the MFIR domain by Spitzer, with both IRAC and MIPS (see Section 3, and Appendix B for details). 

By means of the M R - M B h relation, Wu (2009b) estimated a central BH mass of 5.012 x 1O 8 M . 

In the X-ray domain, the source was observed by ASCA and subsequently by XMM-Newton (Guainazzi et al. 2004, 
and references therein): it displayed the typical features of obscured AGN, despite the optical classification as a 
Seyfert 1; the detection of a prominent iron K a emission line strongly suggested a Compton-rcflcction scenario for the 
X-ray emission, but the spectrum could be satisfactorily fitted by several different models, and a contribution from 
upscattered IR radiation off the lobes could not be ruled out by the authors. Depending on the assumed model for the 
X-ray spectrum, the luminosity of the source in the 2-10 keV range varies from 4 x 10 42 erg s _1 to > 9 x 10 43 erg s _1 . 

References for spectrum and SED of Fig. 1 

Tinti et al. (2005); Dallacasa et al. (2000); Knapp et al. (1990); Stanghellini et al. (1993); Guainazzi et al. (2004); 
this work (Spitzer/IRAC, and Spitzer/MlPS data); NED. 

IERS B2128+048 (PKS 2127+04; OX +046) 

This is the most radio powerful source of our sample, with P 5GHz = 10 27 - 8 W Hz" 1 . Early VLBI observations at 18 
cm were made by Hodges et al. (1984): they found the radio source composed of two resolved components separated by 
29 mas (225 pc). More recent VLBI mapping at 6 cm by Stanghellini et al. (1997) revealed a radio emission extending 
for ~35 mas in NW-SE direction, with the typical scaled-down structure of a Classical Double. At the edges, there are 
two hotspots/lobes; the central, weak component is assumed to be the core; a jet connects the core to the northern 
lobe. At 15 GHz, the VLBA image shows an aligned structure with four main components, the outer of which are 
more resolved than the inner ones, although the southernmost is rather weak (Stanghellini et al. 2001): based on 
morphology, the source is classified as a CSO, with the hotspots being associated with the outer components. 

The source integrated radio spectrum is convex, fairly stable, and GHz-peaked, and the classification of the object 
as a GPS source (O'Dea 1998, and references therein) was confirmed also recently by Torniainen et al. (2007). 

The optical counterpart of this radio source is a very faint (r = 23.3) and red (r — i = 1.85) galaxy (Biretta et al. 
1985; de Vries et al. 2000), located at z = 0.99 (Stickel et al. 1994). The redness of the spectrum is confirmed by de 
Vries et al. (2000), who also reported the remarkable brightness in the NIR K band. 
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By means of the Mr — Mbh relation, Wu (2009b) estimated a mass of 1.698 x 1O 8 M for the central BH. 

In the X-ray regime, the source was pointed at by Einstein in 1979 (Ledden & O'Dcll 1983), but the low photon counts 
only enabled the estimate of an upper limit to its flux in the 0.23-3.77 keV band. The source was then successfully 
observed by Chandra in 2002 (Guainazzi et al. 2006; Siemiginowska et al. 2008) with a 2-10 keV luminosity of 4.4 x 10 44 
erg s~ 4 ; the data spectral fit gave a slope a = 0.5^^7 and a column density A r H=3.0l 8 ;j! ) x 10 21 cm~ 2 . 

References for spectrum, and SED of Fig. 1 
de Vries et al. (2000); Snellen et al. (1996); Guainazzi et al. (2006); Siemiginowska et al. (2008); NED. 

IERS B2352+495 (COINS J2355+4950; OZ +488) 

This radio galaxy has a radio power Psghz = 10 26 3 W Hz -1 . It was imaged at 4.99 GHz by Pearson & Readhead 
(1988) and Conway et al. (1992), who reported on a complex structure surrounding a compact core: two distinct 
components were found, one with a simple structure, and the other with a complex structure. Imaging at 1.67 GHz by 
Wilkinson et al. (1994) also showed a two-sided morphology, which led the authors to classify the source as a CSO. No 
radio structure on scales larger than 0.2" (7 kpc) has been detected for this source (Perley 1982). Consistently with 
Pearson & Readhead (1988) and Conway et al. (1992), VLBI imaging at 4.99 GHz by Fey et al. (1996) confirmed the 
complexity of the source structure, revealing multiple components dominated by a compact core, a region of ~2 mas 
(<~7 pc); two components within the central core, with flux ratio 2:1, contain ^75% of the total flux density at 4.99 
GHz. 

The source exhibits weak variability (15%) on time scales of a few months (Seielstad et al. 1983) and 0.2%-0.7% 
polarization at 5 GHz (Perley 1982) The radio spectrum is convex and GHz-peaked, and the object's classification as 
a GPS source (O'Dea 1998) was confirmed also by Torniainen et al. (2007). 

The source is hosted by a galaxy with redshift z — 0.237 (O'Dea et al. 1991). In the NIR band, the source is one of 
the few displaying a higher nucleation (lower extended emission) in the J band than in the H and K bands (de Vries 
et al. 1998a). 

Wu (2009a) and Wu (2009b) recently estimated a central BH mass of 1.58 x 1O 8 M and 1.698 x 10 8 M Q by means 
of the Mbh — o~ and Mr — Mbh, respectively. 

IERS B2352+495 was firstly observed in the X-rays by ROSAT and ASCA, being detected by neither of them (O'Dea 
et al. 2000). Observations by XMM-Newton in 2004 did reveal the source as a 2 — 10 keV emitter with a moderate 
luminosity of 4.6 x 10 42 erg s _1 (Vink et al. 2006) . The photon statistics was not good enough to independently 
estimate spectral index a and column density Ar . By fixing a — 0.75, Ar= (0.66 ± 0.27) x 10 22 cm -2 was obtained. 

A search for 21-cm absorption successfully detected two features, yielding Ari of 0.28xl0 20 cm~ 2 and 
2.56xl0 20 cm~ 2 (Vermculcn et al. 2003), assuming T s = 100 K and unitary covering factor. Recently, Araya et al. 
(2010) performed high angular resolution Hi measurements of the source, detected the two above absorption features 
against the radio jet, and associated the broader of them with a circumnuclear gas cloud. 

References for spectrum, and SED of Fig. 1 
de Vries et al. (1998a); Snellen et al. (1996); Vink et al. (2006); NED. 

ANALYSIS OF THE SPITZER DATA 
Analysis of the IRAC data 

We analysed the IRAC archival data of IERS B1345+125 and IERS B1404+286 in the following way. The pipeline 
calibrated BCD files were recombined with MOPEX rebinned to l/5th the native pixel size of 1.2" in an attempt 
to separate out contribution from surrounding diffuse emission seen in higher resolution NIR images. The smallest 
practical aperture (r = 2.4") was used to extract fluxes most indicative of the central point source, with aperture 
corrections applied. An uncertainty of roughly 3% in the absolute calibration (Reach et al. 2005) was added in 
quadrature with measurements of the uncertainties due to background fluctuations and the Poisson noise. 

Special care was taken for a proper background evaluation. In the case of IERS B1404+286, in order to avoid the 
diffraction spike from a nearby point source ~83" to the SW of the GPS source, half of an r = 14.4" — 24" annulus was 
used to measure the background. In the case of IERS B1345+125, due to obvious contamination from field sources, 
the background was estimated from circular apertures in source-free regions. 

To quantify contributions from extended emission surrounding the point source in the IRAC images, corresponding 
fluxes in a larger (r = 12") aperture were measured. In IERS B1345+125, this formally resulted in 30%, 16%, 6%, 
and 6% larger fluxes than those in the original r = 2.4" aperture at 3.6, 4.5, 5.8, and 8.0 /xm, respectively. This 
confirmed the presence of diffuse emission surrounding the point source (already noticeable by eye) in the two shortest 
wavelength (highest-resolution) bands, but we consider any similar evidence in the longer wavelength images to be 
marginal at best due to the large aperture corrections (37% — 57%) applied to the small aperture fluxes. 

A similar test of IERS B1404+286 hints at an additional (5%) extended flux in the 3.6 /im image, but the contribution 
is < 2 — 3% in the other images. 

Analysis of the MIPS data 

We analysed the MIPS archival data of IERS B1345+125, IERS B1404+286, and IERS B0500+019 as follows. 
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Using a standard r = 35" (r = 30") circular aperture and a background annulus of 40"-50" (40"-60") for the 24 
/zm (70 /Ltm) image, we measured aperture-corrected source fluxes on the pipeline calibrated PBCD mosaic images. 
The flux conversion utilized the latest (S18) 17 conversion factors of 0.0454 and 702 MJy/sr per instrument unit, with 
absolute uncertainties of 4% and 7%, respectively. As in the IRAC data, these systematic uncertainties were added 
in quadrature with the background and counting noise. Additionally, the aforementioned point source 83" from IERS 
B1404+286 is unresolved from the GPS source in the lower resolution IRAS images, so presumably contaminates the 
older measurements (we measured 117 ± 5 mJy at 24 /im). 
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[10 erg s J 


IERS B0026+346 


0.517 


35 


211.4 


6.04 


27.1« 


630.96 1 


957.17 






1.84 2 


IERS B0108+388 


0.66847 


6 3 


40.9 


6.82 


27.2 3 


3900-4760 4 ' 5;6 


6510-8340 


0.25±0.04 7 


367 ± 48 s (k) 
310 ± 70 7 (k) 


2.00 9 






15 3 














417 10 (k) 




IERS B0500+019 


0.58457 


96.3 


6.42 


27. 1 3 


2000 4 


3169.14 






4.95 2 


IERS B0710+439 


0.518 


25 3 


151.1 


6.04 


27.0 3 


1500-1900 4 ' 5 


2277-2884 


0.37±0.17 7 


550 ± 160 7 (k) 
932 10 (k) 


2.15 9 




















1200-1800 n (s) 




PKS B0941-080 


0.228 


50 3 


177.4 


3.55 


26.0 3 


500 4 - 5 


614 




0.00578 2 , 0.005 12 
0.33 9 


IERS B1031+567 


0.450 


40 3 


224.0 


5.60 


26.8 3 


1300 - 1400 5 ' 4 


1890-2030 


0.68±0.15 7 


620 ± 140 7 (k) 
1836 10 (k) 


IERS B1345+125 


0.12174 


80 h ' 3 


170.4 


2.13 


25.9 3 




673 






0.40 14 , 0.41 12 


IVS B1358+624 


0.4310 


70 3 


382.3 


5.46 


26. 8 3 


500 4 ' 5 


715 






2.538 9 


IERS B1404+286 


0.07658 


7 3 


9.9 


1.41 


25.4 3 


4900-5340 15 ' 6 


5280-5750 


0.28±0.11 16 


320±210 17 (k) 


0.038 - >0.9 18 




















224 1(, (k) 


4.20 2 , 3.21 12 


IERS B2128+048 


0.990 


35 1 


271.9 


7.77 


27.8 3 


600-800 4 * 5 


1194-1592 






IERS B2352+495 


0.23790 


50 3 


183.2 


3.66 


26.2 3 


700-900 5 ' 4 


870-1110 


0.40±0.13 7 


1200±400 7 (k) 


0.0488 9 



3003 10 (k) 



References. — (1) Snellen ct al. (1996); (2) Guainazzi ct al. (2006); (3) O'Dea (1998); (4) dc Vrics et al. (1997); (5) Stanghcllini ct al. (1998); (6) Tinti ct al. (2005); (7) Taylor 
ct al. (2000); (8) Owsianik ct al. (1998); (9) Vink ct al. (2006); (10) Polatidis & Conway (2003); (11) Rcadhcad ct al. (1996); (12) Siemiginowska ct al. (2008); (13) Stanghellini 
ct al. (2001); (14) O'Dea ct al. (2000); (15) Dallacasa ct al. (2000); (16) Stanghcllini ct al. (2002); (17) Liu ct al. (2000); (18) Guainazzi ct al. (2004). 

a Redshift, taken from NED. 

^Overall angular size of the radio source, unless otherwise stated. 
c Observed turnover frequency (or frequency range) of the radio spectrum 
^Intrinsic (source-frame) turnover frequency (or frequency range) of the radio spectrum 
c Hotspot separation velocity, derived from the projected angular velocities. 
^Source age, estimated by means of kinematic (k) and/or spectral ageing (s) methods. 
s Dcrivcd from Fsghz- 
The size refers to the one-side, southern radio jet. 
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TABLE B2 

Flux densities of the sample's sources observed by the Spitzer Space Telescope 



Source name 




IRAC 






MIPS 




-F3.6 (im 


-F4.5 fJ,m -F5.8 fim 




-F24 pm 






(mjy) 


(mjy) (mjy) 


(mjy) 


(mjy) 


(mjy) 



IERS B0500+019 ■ • ■ ■ • ■ • ■ • • • ■ 6.3 ± 0.3 

IERS B1345+125 8.2 ± 0.6 13.4 ±0.9 20.9 ±1.1 42.3 ±1.9 479 ± 21 1947 ± 186 

IERS B1404+286 37.0 ± 1.7 47.0 ± 2.0 54.5 ± 2.3 78.7 ±3.1 382 ± 17 804 ± 102 



Note. — 



Details on the analysis of the Spitzer data can be found in Appendix B. 



TABLE B3 

SED MODELING: Fit parameters and derived quantities 



Source name Lj L uv L m L opt si s 2 i/ p u m Lj/(10L UV ) ^uv/^Edd [^BH,s] (a) 

(10 45 erg s- 1 ) (10 45 erg s" 1 ) (10 45 erg s" 1 ) (10 45 erg s" 1 ) (GHz) (THz) 



IERS B0026+346 


0.70 


10 


6.00 


1.00 


0.70 


2.0 


0.7 


5 


0.007 


0.390 


1.698*] 


IERS B0108+388 


2.05 


5 


0.50 


0.60 


1.80 


3.2 


6.0 


5 


0.041 


0.417 


0.794*] 
1.698*] 


IERS B0500+019 


1.80 


5 


3.50 


1.00 


1.70 


2.0 


2.0 


5 


0.036 


0.195 


IERS B0710+439 


2.80 


10 


4.30 


0.90 


2.00 


2.1 


2.5 


5 


0.028 


0.264 


2.512*] 


PKS B0941-080 


0.40 


1 


0.08 


0.80 


2.00 


2.5 


0.5 


5 


0.040 


0.132 


0.501*] 


IERS B1031+567 


0.87 


5 


1.40 


0.30 


1.90 


2.0 


1.0 


5 


0.017 


0.263 


1.259*] 


IERS B1345+125 


0.17 


5 


6.00 


0.90 


1.40 


2.0 


0.4 


5 


0.003 


0.506 


0.654*] 


IVS B1358+624 


2.10 


5 


6.00 


0.35 


1.75 


2.7 


0.6 


5 


0.042 


0.209 


1.585*] 


IERS B1404+286 


0.30 


2 


0.70 


0.50 


2.15 


5.6 


5.0 


13 


0.015 


0.026 


5.012*] 
1.698*] 


IERS B2128+048 


4.00 


5 


3.00 


0.65 


1.70 


2.5 


1.2 


5 


0.080 


0.195 


IERS B2352+495 


0.45 


5 


0.60 


0.35 


2.00 


2.0 


1.0 


5 


0.009 


0.132 


2.512*] 



Note. — The fit parameters (cols. 2—9) are described in Sections 2 and 4; the derived quantities (cols. 10 and 11) arc discussed in Section 6. Parameters common to all 
the SED fits, not included in this Table, arc as follows: n = 0.1 cm" 3 , 7 min = 1, 7 max = 100 m p /m C) r? e = 3, r/ B = 0.3, z> opt = 2.0 x 10 14 Hz, y uv = 2.45 x 10 15 Hz. We 
assumed hotspot propagation velocity v-^ — 0.1c when no hotspot separation velocity Uh,sep was available (see Tabic 1); from the known i^sep; we derived t>h by using the 
appropriate composition law for velocities, yielding in most cases C=L 0.5^h,scp.- 

( a )The Eddington luminosity was computed according to the relation: L-^dd — 1.51 X 10 46 Mbh/(10 8 Mq) erg s~ 1 (e.g., Krolik 1998). The assumed BH mass is given in square 
brackets in units of 10 s Mq: it is the most recent estimate of the source BH mass when available (from Wu (2009a) when marked with from Wu (2009b) when marked 
with and from Dasyra ct al. (2006) when marked with °), and the average value of the BH mass estimated by Wu (2009b) for his GPS-source sample (Mbh,8 — 1-698) 
otherwise (i.e., for masses marked with *). 
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TABLE B4 

X-RAY SPECTRAL PARAMETERS AND RADIO Hi COLUMN DENSITIES. 



Source name 






r 


Ref. 




Ref. 




(10 20 cm- 2 ) 


(10 22 cm- 2 ) 






(10 18 x (T e /c f ) cm- 2 ) 




IERS B0026+346 


5.6 


1 0+ ' 5 
1 - u -0.4 


-, 4 o+0.20 
1 -^°-0.19 


2 






IERS B0108+388 


5.8 


57 ± 20 


1.75( c > 


3 


80.7 


8 


IERS B0500+019 


8.3 


u.a_ 2 


1 62+ ' 21 
1 - DZ -0.19 


2 


6.2 (4.5,2.5) 


8 


IERS B0710+439 


8.11 


0.44 ± 0.08 


1.59 ± 0.06 


3 






PKS B0941-080 


3.7 




2( c > 


2 


<0.80( d > , <1.26( c ) 


6, 10 




3.7 




z - uz -1.03 


9 








3.67 


< 1.26^ e ^ 


2 - 28 -0.61 


9b 








3.67 


< 0.53 (o) 


(1.7 - 1.9) (°) 


9b 








n e:ft 


n ten i nig 




Q 
O 


<U.O y , <±.^D^ ' 


ft in 
0, 1U 




1.1 


4 2+ 40 


1 6+ 1 ' 2 
i - D -0.8 




6 . 2 ( f ) 


7 




i.y 


„ -,,+0.636 


-1 1n +0.29 


y 






IVS B1358+624 


1.96 


3.0 ± 0.7 


1.24 ± 0.17 


3 


1.88^ 


6 


IERS B1404+286( h ) 


1.4 


n ,0+0.12(1) 


2.lt°; 6 3 « .7l° i«) 


4 


1.83( g > 


6 






0.12t°;0|W , 24.01^ °(J) 


„ „_i_n 4 












o.ii ±o.o5( i - k ) 


2 21 +0.19(1) 












io+0 l3(i) 


2.6±0.5W, 20' c ) 












u - uy -0.06 


2.2 ±0.4 








IERS B2128+048 


5.2 


n o+o.8i 


i r+0.6 


2 








5.2 




1 28+ ' 42 


9 






IERS B2352+495 


12.4 


0.66 ±0.27 


1.75< c > 


3 


0.28; 2.56<S' n ) 


6 



References. — (1) O'Dea et al. (2000); (2) Guainazzi ct al. (2006); (3) Vink ct al. (2006); (4) Guainazzi ct al. (2004); (6) Vermeulen 
ct al. (2003); (7) Mirabel (1989); (8) Carilli ct al. (1998); (9) Siemiginowska ct al. (2008); (9b) Siemiginowska 2008, priv. comm.; (10) 
Pihlstrom ct al. (2003). 



( a ) Equivalent total hydrogen column density, assumed to be located at the source rcdshift; the values in boldface are those used for the 
correlation analysis (see Section 5.2, and Tables B5 and B6). 

(k) Neutral hydrogen column density derived from low resolution measurements, unable to determine the location of the absorbing gas; 
the values in boldface are those used for the correlation analysis (see Section 5.2, and Tables B5 and B6). 
( c ) Fixed. 

^Derived from the value of A^hi,2ct given by (6), which was obtained from the relation: 1.82 X 10 18 X T s X T2 CT X AV, under the assumption 
T s — 100 K, AV — 100 km s _1 , and cj — 1; no detection of absorption lines. 
^ c ^3cr upper limit. 

( f )iVHi~ 10 22 x (T s /100 K) cm 2 was obtained by Morganti et al. (2004b) from high-resolution measurements that constrain the location 
of the absorbing gas. 

Derived from the value of -/V H i,2(t given by (6), which was obtained from the relation 1.82 X 10 1S xT s x r 2cr X AV, under the assumption 
of T B = 100 K, and c f = 1. 

( h -*For this source, we report the X-ray spectral parameters derived by Guainazzi ct al. (2004) by means of different models for the X-ray 
emission; the first three parameter sets refer to the bow-tics represented in Fig. 1; the third set is bes fitted by our model. 
WSoft X-rays (E > E hlo&k ). 
(j)Hard X-rays (E > £ brcak ). 

{ -'Associated with a Compton-rcflcction model for the hard X-rays, with Af H hard > 9 x 10 23 cm 
(^Intrinsic, for the Compton-reflection model. 
( m )3tr errors. 

'■"•'After submission of this paper, Araya et al. (2010) published high resolution measurements that constrain the location of the absorbing 
gas, and derived the values N U i= 7.3 X 10 21 cm -2 and AT H i = 5.2 X 10 22 cm -2 , by assuming T s = 8000 K. 



X-RAY GPS/CSO GALAXIES: MODELING THE SEDs 



TABLE B5 

Nfi-Nfil CONNECTION: CORRELATION AND REGRESSION ANALYSIS FOR DETECTIONS 



Sample 3 


N 


Pearson 


Spearman 


Kendall 


Linear reg 


;rcssion 






r Prob. b 


p Prob. b 


t Prob. b 


Slope 


Intercept 


D5 


6 


0.9966 1.739c-05 


0.3947 0.4387 


0.2981 0.4008 


0.837±0.109 


5.151±2.269 


D5+AD 


7 


0.9967 1.229e-06 


0.6301 0.1294 


0.5143 0.1048 


0.981±0.100 


2.122±2.067 



a D5 is our sub-sample of Nh-Nhi detections; AD indicates Nh-Nhi detections for additional sources reported by 
Tengstrand et al. (2009); sec Section 5.2 and Table B4 for details. 

^Probability of the null hypothesis of no correlation being true. It is a two-sided significance level: because we 
are looking a priori for a positive correlation, this value should actually be divided by 2, improving the significance 
by a factor of 2. 



TABLE B6 

Nji-Njii CONNECTION: CORRELATION AND REGRESSION ANALYSIS FOR DETECTIONS AND UPPER. LIMITS 



Sample a 


N 


Generalized Spearman 


Generalized Kendall 


Schmitt's linear regression 






p Prob. b ' c 


S d Prob. b 


Slope Intercept 


D5+U2 


8 


0.670 0.0764 


26. (2=1.749) 0.0804 


0.9433±0.0621 2.8378±1.3403 


D.5+U2+ADU 


11 


0.738 0.0197 


54. (2=2.402) 0.0163 


1.0711±0.0034 0.0775±0.0695 



a D5 is our sub-sample of N^-Nm detections; U2 is our sub-sample of JVh and/or -/Vhi upper limits; ADU is the set 
of detections of and upper limits to, TVh and/or -/Vhi for additional sources reported by Tengstrand ct al. (2009); see 
Section 5.2 and Table B4 for details. 

b Probability of the null hypothesis of no correlation being true. It is a two-sided significance level: because we are 
looking a priori for a positive correlation, this value should actually be divided by 2, improving the significance by a 
factor of 2. 

c According to the ASURV Rev. 1.2 software package, this value is accurate only if N > 30. 
d S is denned in Isobc ct al. (1986). 



